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ABSTRACT 


ObSGirvatioii of Bow Shock Protons at the Lunar Orbit 

by 

John Louis Benson 

Protons with energies ranging from about 500 eV to 
3500 eV are observed by the Suprathermal Ion Detector Ex- 
periment on both the dusk and dawn sides of the magneto- 
sphere . On each lunation these particles appear as a rather 
continuous phenomenon for 3 to 5 days after crossing from 
the dawn-side magnetosheath into the solar wind and for 
about 2 days prior to entering the dusk-side magnetosheath. 

Transverse proton flows of this type have been observed 
by AiSbridge et al.^ [1968] , Scarf et al . c [1970] and are be- 
lieved to be deviated solar wind particles that are acceler- 
ated at the earth’s bow shock and emitted outwardly along 
interplanetary field lines at greater than solar wind 
energies . 

Data from the SIDE and from the Explorer 35 lunar orbit- 
ing magnetometer have been analyzed and these data indicate 
that the transverse ion flows observed by the SIDE in the pre 
and post bow shock crossing regions of the lunar orbit are due 
to these deviated solar wind particles. 

A computer model based on EXB drift trajectories for 
particles leaving the shock has been developed and synthetic 
particle data produced by this model are in good agreement 
with the observed data. 
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Chapter 1 

This thesis is concerned with the study of the plasma 
environment of the moon; particularly when the moon is outside 
the earth's magnetosphere and in the interplanetary medium. 
Much has been learned about this environment as a result of 
placing particles and fields experiments directly on the lun- 
ar surface during the six Apollo moon landings and the follow- 
ing work consists of the analysis of data returned to earth 
by one such experiment, the Suprathermal ^on Detector Experi- 
ment, otherwise known as SIDE. 

The plasma environment of the moon should not be consid- 
ered independently of the earth's magnetosphere and the in- 
terplanetary medium and a brief review of the basic ideas and 
observations concerning the complex interactions of the solar 
corpuscular radiation, (solar wind) , and the near extra-mag- 
netospheric region with the earth moon system is in order 
before going on to a detailed discussion of the results of 
this analysis. 

1.1 General Survey 

The s\in continuously ejects a tenuous supersonic plasma, 
known as the solar wind, that consists mostly of protons and 
electrons in equal amounts . The plasma is tenuous in that 
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typical number densities are on the order of a few particles 
per cubic centimeter and supersonic in that the bulk flow 
velocity of the plasma is greater than the Alfv^ speed in 
the medium. The Alfven speed is, 

=B/C4 7t p ~ 50 km/sec 1-1 

in the interplanetary mediim and is the quantity analogous to 
the sound speed in a collision dominated medium. The average 
bulk flow velocity of the solar wind is about 400 km/sec 
giving an Alfvenic Mach number, about 8, indicating 

the highly supersonic flow. 

The solar wind drags the sun's magnetic field lines 
along with it due to the essentially infinite conductivity 
in the plasma in such a manner that the ambient interplaneta 3 :y 
magnetic field lines intersect the earth's orbit at approxi- 
mately a 45° angle to the earth-sun line. The field lines 
in figure 1-1 are shown at this "garden 'hose angle." 

The radial flow of plasma from the sun appears to come 
from a direction about 5° to the west of the solar direction 
due to the earth's orbital velocity of about 30 km/sec. 

Upon reaching the orbit of the earth, the solar wind 
plasma is unable to penetrate the earth's dipolar field; it 
therefore compresses the field creating a cavity in the flow 
known as the magnetosphere. The solar wind flows around this 


Figure I-l. View of the magnetosphere and the 
lunar orbit in the ecliptic plane. Typical magnet- 
ic field lines in the interplanetary medium and in 
the magnetosphere are shown. 
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cavity in a fashion resembling hydrodynamic flow. The 
boundary of the cavity is called the magnetopause and the 
distance from the center of the earth along the earth-siin 
line to this boundary is on the order of 10 (earth radii) , 

The supersonic nature of the solar wind gives rise to a 
detached shock wave a fev/ R^ upstream of the magnetopause . 

This standing shock/ known as the bow shock, v;as predicted by 
Zhigulev and Romishevskii in 1959 and by Axford in 1962. Its 
existence was confirmed by the IMP 1 satellite, Hess [1964] . 
Due to the excessively long collision mean free path of a 
solar wind particle, — 1 A.U., the shock is not collision 
dominated or even collision influenced. Nevertheless, the 
high conductivity in the plasma causes it to behave somewhat 
like a fluid and shock calculations using the "fluid dynamic 
analogy", Spreiter et al . [1966], work fairly well in pre- 
dicting the location of the shock. See figure 1-2. A few 
R^ upstream of the shock there is only relatively cool, 

(T ~ 10 °K) , tenuous, (n 10/cm ) , unshocked plasma with 

a weak interplanetary field, (B ~ 57 ), imbedded in it. 

For a stimmary of the properties of the quiet solar wind, see 
table 1-1. Closer to the shock, where interplanetary field 
lines intersect the shock surface, complex interactions of 
the solar wind and the shock take place and energetic parti- 
cles have been observed flowing upstream of the shock in this 
region by several investigators; Asbridge et al. [1968], 
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Figure 1~2 . Measurements of the positions of the mag- 
netopause and bow shock frora IMP 1, [after Ness^ (1965)]. 



Figure 1-2 
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Table 1-1 

Observed Properties of the Quiet Solar Wind 


bulk flow velocity 
number density 
proton temperature 
proton thermal anisotropy ratio 
electron temperature 
magnetic field strength 
proton flux 
kinetic energy flux 
kinetic energy density 
proton thermal energy density 
electron thermal energy density 
magnetic field energy density 


320 km/sec 
8 cm ^ 

4X10'^°K 

2 

1 to 1.5X10^ ‘’K 
57 


2.4X10® cm ^sec"^ 
-2 “1 

.22 erg cm sec 
7X10 ^ erg cm ® 
6X10"^^ erg cm'^ 

1.5X10 erg cm 

T^-10 -3 

10 erg cm 


3 
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Prank and Shope I1968] , Anderson [1968]and Freeman I1972J . 

I 

These energetic particles, specifically protons in the energy 

range of .5 to 3.5 keV are the siibject of this thesis. A 

more detailed discussion of these observations is given belov/. 

After passing through the shock, the temperature of the 

plasma is increased, {T ~10^“K) , the magnetic field is com- 

-3 

pressed, (B ~ 20 y ) , density goes up, (n ~ 30/(.nn ) and the 

bulk flow velocity goes down. Some of the ordered streaming 
energy is converted to disordered thermal energy . Figure 1-3 
summarizes the changes in plasma parameters across the shock. 
The region of space between the magnetopause and the bow 
shock is known as the transition region or the magnetosheath 
and the detailed plasma properties of this region near the 
lunar orbit are discussed by Fenner 11971,1974]. 

The inside of the magnetosphere is essentially shielded 
from the solar wind flow and figure 1-4 depicts the general 
particle and field morphology relevant to the Lxanar environ- 
ment in this region of space. The earth’s dipolar field is 
compressed on the day side in to about 10 and drawn out on 
the night side at least several tens of past the lunar 
orbit iNess et al. , 1967] , in the anti— solar direction. 

The dipolar field lines are closed near the earth, but 
the distance in the anti-solar direction to which this clos- 
ing extends is not well established. At a distance of more 

than about 10 R , the field tends to point either toward or 
e 
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Figure l-3a. Typical changes in plasma parameters 
across the bow shock. Note that C-soiand speed and 
b=Alfven speed. =mach number =Alfven mach number, 

[after McKenzie, (1970) ] . 

Figure l-3b. Measured changes in plasma parameters 
across the shock, [after Wolfe and McKibbin, (1968)]. 

Figure 1-4. View of magnetosphere in noon-midnight 
meridian plane, [after Ness, (1965)]. 
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away from the sun. Although the field lines may be v^eakly 

closed CB„ << 1) at 6Q R , the morphology of the geomag- 

ss 

netic tail in this region can generally be described as a 

planar neutral sheet v/ith the field lines roughly parallel to 

the sheet above and below v/ith the sense of the field being 

toward the sun north of the neutral sheet and away from the 

sun south of the neutral sheet . The thickness of the neutral 

sheet at 60 R is not well known, but Bowling, [1973] has set 
e 

a value of about 3 R^ at a distance of 30 R^ in the anti- 
solar direction. 

The magnetic neutral sheet region is also filled with a 

3 

relatively warm, (T^ ~ 10®°K) , rarefied, {n - 10 /cm ) , 
plasma that surrounds the neutral plane with a thickness of 
about 3 to 7 R at the lunar distance as determined by CPLEE , 
[Rich, 1973] . See chapter 2 for a description of the Charged 
Particle Lunar Environment Experiment. The source of the par- 
ticles in this "plasma sheet" is not well established, but it 
is believed that they are solar wind particles that have 
entered the magnetosphere through the magnetopause or through 
open tail field lines at large distances in the anti-solar 
direction. 

The above brief description of the regions of the magnet- 
osphere and interplanetary medi\m through which the moon pass- 
es allows one to summarize the changing plasma environment of 
the moon , disregarding any lunar atmospheric phenomena , as 


follows . 
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Beginning at new moon, Csee figure 1-1) , the moon is in 
the undisturbed solar wind where the ordered streaming plasma 
rams the stinward side directly without going through a shock 
wave and the earthward side faces a plasma void that lines up 
with the solar wind aberration direction. The absence of a 
lunar bow shock was discovered from Explorer 35 magnetometer 
data, [Colburn, 1967] and the plasma void was identified using 
a plasma probe flown on the same spacecraft , [Lyon et al . , 

1967] . The Explorer 35 magnetometer has also determined that 
a slight enhancement of the interplanetary field occurs while 
the satellite is in the moon's optical shadow, [Colburn et al., 

1968] . 

At first quarter the conditions are essentially the same 
as at new moon, as long as the interplanetary field does not 
change its streaming angle appreciably. That is, at first 
quarter the interplanetary magnetic field lines still do not 
connect the moon and the bow shock so that the bow shock can 
not yet influence the lunar environment . 

When the moon proceeds toward the dusk magnetosheath, it 
will begin to encounter field lines connected to the bow shock 
with the probability of observing bow shock associated part- 
icles becoming larger as the moon moves closer to the shock. 

After passing into the dusk magneto sheath, the plasma 
environment changes considerably from that just outside the 
shock. As discussed above, the density increases, the bulk 
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flow velocity decreases and the temperature increases. The 
magnetic field is enhanced in this region and particle spec- 
tra become very stable with peak energies aroiind 500 eV. 

By the time the moon has reached full phase it is well 
into the geomagnetic tail where it experiences multiple 
.neutral sheet crossings and encounters v^ith the plasma sheet, 
densities decrease from raagnetosheath values and temperatures 
increase. 

As the moon passes from the full to the last quarter 
phase it crosses the dawn magnetosheath whexa the plasma par- 
ameters are similar to those in the dusk magnetosheath but 
differences in turbulent effects are observed. The dusk 
magnetosheath exhibits a correlation in turbulence and geo- 
magnetic activity (Kp) and the dawn magnetosheath does not. 

[Fenner, 1974] . 

The region between the dawn bow shock and the point at 
which the last interplanetary field line tangent to the shock 
intersects the lunar orbit, (see figure 1-1), is the primary 
region of interest in this thesis since the probability of 
observing bow shock protons is greatest here. 

1.2 Observations of bow shock associated phenomena, upstream 
of the shock . 

Several investigators have observed energetic particles 
and wave disturbances upstreani of the earth's bow shock. A 
summary of those observations follows. 
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Waves of frequency 10-100 Hz propagating along 
interplanetary field lines have been detected upstrea of the 
earth's bow shock by the OGO-1 search coil magnetometer, 

[Smith et al,, 1967]. These waves are believed to be precur- 
sor waves propagating upstream in the whistler mode parallel 
to the field lines. Similar observations of whistler waves 
upstream of the shock were made by Heppner et al. [1967] 
although the frequencies were about 1 Hz with the waves being 
detected in small packets of 4 to 6 cycles . The amplitudes 
of these waves were found to diminish with increasing distance 
from the shock which led Heppner to suggest that they were 
generated at or were associated with the shock itself . 

A comprehensive study of upstream waves was made by 
Fairfield [1969], using Explorer 34 data. He determined that 
waves of frequency 0.01 to 0,05 Hz upstream of the bow shock 
are a permanent feature of the region of space extending 
from the bow shock out to about 34 E,^. These waves are 
observed 90% of the time when the interplanetary field inter- 
sects the shock and Fairfield believes them to be locally 
generated Alfven waves due to the outward streaming protons 
observed by Asbridge and other investigators . He argues that 
the whistler waves must have a fairly high frequency 
( V ^ 10 Hz) to propagate upstream as shown by Smith et al. 
[1967] . Processes that would explain the necessarily large 



13 


Doppler shifting to .01 Hz in the rest frame of the spacecraft 
v 7 ould be difficult to explain. 

Greenstadt et al. [1968] have contributed further to the 
association of these low frequency waves with particles by 
showing a correlation in magnetic field oscillation and 
particle fluxes outside the shock. Energetic electrons 
(energy ^40 keV) have also been observed in the upstream re- 
gion by Anderson [1968], and by Jokipii [1968]. Protons of 
1 keV were observed upstream of the shock by Frank and Shops 
[1968] , with the IMP-4 satellite v?hile the detector was look- 
ing in the anti-solar direction. Freeman [1972] has also 
observed energetic ions in the deep lunar night with the 
suprathermal ion detector on the liinar surface when the detec- 
tor was looking as near as 17° to the anti-solar direction. 
Freeman points out that these ions are probably not from the 
earth's bow shock and that they are most likely associated 
with deviated solar wind flow around the moon. Examples of 
outwardly f levying protons are reported by Asbridge et al . 
[1968], using data from six Vela satellites. The results 
reported in their paper are particularly relevent to the 
present work and this paper will now be discussed in more de- 
tail. 

Asbridge and co-workers report the results of the 
analysis of particle data from the six Vela satellites which 
are in nearly circular earth orbits at 18 R^. The energy 
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analyzers on hoard are capable of detecting protons in the 
range of 0.1 to 20 keV. All of the Vela satellites are spin 
stabilized and sample only certain directions in spacecraft 
azimuth. Vela's 2 and 3 spin at about 2 rps and take five 
measurements near the solar direction and others at 90°, 180°, 
and 270° from the solar direction. The Vela 4 satellites spin 
at about 1 rpm taking seven measurements of proton spectra 
near the solar direction and a measurement each 22 1/2 of 
rotation thereafter. 

Several examples of proton flows in the energy range 3 to 
6 keV are presented in this paper. These protons are observed 
to come from the direction of the earth's bov; shock and figure 
l-5a,b,c shows two examples of the dawn side flows and one 
example of a dusk side flow. An analysis of four months of 
Vela 3A and 3B data shows that these outwardly flowing protons 
are a permanent feature of the dusk-davm bow shock region 
and figure l-5d summarizes these results. Arrows indicate 
direction of flow and dots indicate flow in the solar direc- 
tion. A higher incidence of outwardly flowing protons is 
noted when the spacecraft are on the dawn side .of the magnet- 
osphere. The authors suggest that these protons are solar 
wind particles that are accelerated at the bow shock and 
emitted outwardly from the shock along interplanetary field 
lines. Figure 1-6 siommarizes the model suggested by these 
investigators and it is shown in this thesis that this model 
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Figure l-5a. Dawn-side transverse proton flow from 
the bow shock as observed by Vela 3B; 291° spacecraft 
azimuth points toward the bow shock while on the dawn 
side, (energy ^ 4 keV) . 

Figure l-5b. Dusk-side transverse flov/ from Vela 2B 

observations at 18 R . (energy 3 keV) . 

0 

Figure l-5c . Vela 3A dusk-side observation of transverse 
flow at 89° from the solar direction, (energy 5 keV) . 

Figure l-5d. Summary of four months of Vela 3A and 3B 
data. Arrows indicate direction of the flow; dots indi- 
cate flow in the solar direction. 

Figure 1-6. Diagram of Asbridga model to explain the 
observed transverse proton flows . 


.Figures 1-5 and 1-6 are from [Asbridge et al. (1968)]. 
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is at least partially correct. A more realistic model is de- 
veloped and discussed along with the Asbridge model in chap- 
ter 3. 


1.3 Research Objective s 

The SIDE observes a variety of ion spectra in the region 
of space between the bow shock and the terminator. The dom- 
inant spectrum is one that is depleted in the lower energies 
and peaks toward the higher energies . Other spectra show 
double peaks with one doubly peaked type occuring near the 
bow shock and the other near the terminator. The lack of a 
complete imders tending of the ion data in this region and the 
fact that the dominant spectrum is what one would expect from 
bow shock protons lead to the following research objectives; 

a. ) to determine v/hether or not the moderate to higher 
energy ions observed by SIDE for several days after passing 
out of the magnetosheath and into the solar wind are parti- 
cles from the earth's bow shock. 

b. ) to identify these particles in the data, further 
explaining some of the more complicated spectra observed in 
this region , iLindeman, 1971] . 

c. ) to propose and develop a theoretical model to ex- 
plain how bow shock protons can reach the lunar orbit. 
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Chapter 2 

2 .1 Basic Theory of instriiment Operation 

e basic idea behind the SIDE and most particle detectors 
with energy discrimination capability is simply that the traj- 
ectory of a charged particle is altered upon its entering a 
region in which there exists an electric field. Detectors 
that are based on this principle often employ a set of parall- 
el plates across which a potential difference is applied to 
generate the electrostatic field. A channel electron multi- 
plier is placed at the end of the plates as a counter. If the 
configuration of the plates is such that the electric field is 
applied perpendicular to the trajectory of the incoming part- 
icle then the amount of deviation from that initial trajec- 
tory depends on the energy of the particle and the strength of 
the electric field. 

The very simplest configuration for such a detector would 
be a pair of straight parallel plates of length L that are 
separated by a distance d. Then if a potential V is applied 
across the plates, a particle entering the region between the 
plates on a trajectory that is parallel to the plates and at 
a distance of d/2 above the bottom plate, must have an energy 
that is at least equal to 

£ = q 2-1 

c 

to get through the plates, where q is equal to the charge of 
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the particle . 

This configuration is not very useful since it only pro- 
vides the lower limit on the particle energy, i.e. a particle 
with energy greater than can also make it through the plates. 
The energies of particles that get into the counter can be 
confined to a narrow range if the parallel plates are curved. 

The SIDE employs such a pair of curved plates in both the 
total ion detector (TID) and the mass analyzer detector (MA) . 

The potential, as a function of r, between the curved 
plates is given by, 

<Afrl= 2-2 

In (b/a) 

where b is the radius of curvature of the outer plate and a is 
the radius of curvature of the inner plate (see figure 2-1), 

V is the total potential difference across the plates. 

The electric field is given by, 



so that the field strength is 

V 


E = 


In (b/a) r 


2-4 


The force exerted on a particle of charge q, moving 
through the curved plates is , 


F = q E 
P 


Figure 2 ~ 1 . Schematic representation of parallel 
plate electrostatic energy analyzers, flat and curved. 
Crossed E and B fields (Wien velocity filter) proceed 
the Low Energy Curved Plate Analyzer to form the Mass 
Analyzer. 





20 


and the centrifugal force on the particle of mass m, is, 


F = 
c 


If these forces are balanced, the particle will make it 
through the plates and into the counter. That is, the energy 
per unit charge of a particle that gets through is a function 
of the potential difference across the plates and the plate 

separation . 

C/„ = V 2-7 

2 ln(b/a) 

Equation 2-7 is only approximately correct since equation 
2-2 was derived for infinitely long cylindrically concentric 
plates. The parallel plates of the detector are merely a sec- 
tion of the concentric cylinders and hence, edge effects are 
important in determining the correct equation relating the en- 
ergy and the potential difference across the plates. The geo- 
metry of the plates and these edge effects give rise to a fin- 
ite range of energies that may be passed at a particular val- 
ue of V. The results of the calibration of the Apollo 14 in- 
strument to determine the size of the energy passband are 
shown in figure 2-2. The bandwidth at the nominal center 
energy is approximately 10% of the center energy. 


2 . 2 Mass Analyzer 


The curved plates of the mass analyzer are preceeded by 


Figure 2-2. Some calibration results of the Apollo 14 
Total Ion Detector. Upper part of the graph shows the 
unidirectional geometric factor as a function of the 
nominal center energy. Lower part shows the energy band 
pass (FWHM) as a function of the nominal center energy. 
Note that the bandpass is approximately 10% of the center 


energy . 


BAND PASS (FWHM) 



Figure 2-2 
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Wien velocity filter which consists of crossed electric and 
magnetic fields that select the incoming particle’s velocity 
before it is energy-analyzed in the curved plates . If the par- 
ticle’s velocity and its energy per unit charge are known, then 
the mass per unit charge is simply, 

ra/q = 2 ^ /{v^q) 2-8 

Upon entering crossed electric and magnetic fields, a 
charged particle experiences two forces, one associated with 
the E field and one associated with the B field. The relative 
strengths of these two fields determines the resultant trajec- 
tory of the particle. The force due to E is. 



and that due to B is, 

\ = qv xt 2-10 

so that if, 

“Pg + = 0 2-11 

the particle goes through the fields undeviated from its ini- 
tial trajectory. This requires that the particle have a velo- 
city, 

v = E/B 2-12 

to pass through the velocity filter and into the curved plate 
analyzer. 
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The E field of the velocity filter is provided by flat 
Parallel plates and the B field is provided by a permanent mag- 
net. By adjusting the voltage on the plates, various ratios 
of *E to^ can be obtained, thus allov/ing only particles of cer- 
tain velocities through the filter. 

The basic idea behind the mass analyzer is also depicted 

in figure 2-1. 

2 . 3 Description of .the instrument 

Figure 2-3 shows schematically the arrangement of the var- 
ious components of the instrument. The SIDE has two positive 
ion detectors, one of which is proceeded by a Wien velocity 
filter, as discussed above, and this detector is known as the 
mass analyzer detector (MA) . The curved plate analyzer, asso- 
ciated with the mass analyzer is connected to a voltage step- 
ping supply and is sensitive to energies per unit charge rang- 
ing from 48.6 to 0,2 eV/q and is referred to as the Low Energy 
curved Plate Analyzer (LECPA) . After passing through the 

velocity filter and the LECPA, an ion enters a counter known 

... 

as a Bendix Channeltron . This channel electron multiplier is 
held at a potential of -3.5 kV to accelerate the ions after 
they pass out of the LECPA. This post acceleration of the ions 
affords a higher detection efficiency of the channeltron which 
emits a cascade of secondary electrons whenever an ion strikes 
its interior. 




■7,'' 


7T 


1 


Figure 2-3 . Schematic diagram of the electrical con- 
nections of the TID and MA. [Hills et al. (1971) ] . 
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The curved plate analyzer associated with the total ion 
detector (TID) , is known as the lU.gh Energy ^urved plate Analy- 
zer (HECPA) . It is also connected to a voltage stepping sup- 
ply and is designed to pass energies per unit charge ranging 
from 3.5 keV/q to 10 eV/q. Ions leaving the HECPA are also 
post accelerated into a channeltron. The physical arrangement 
of the various components of the instrument is shown in figure 
2-4. 

The velocity filter for each of the three instruments 
is different and hence each measures a different mass range. 
Table 2-1 summarizes the energy and mass sensitivities of the 
three instruments and figure 2-5 shows some results of the 
mass analyzer calibration. 

Each SIDE is also equipped with a spider web-like screen 
Csee figures 2 -6, 2 -7) made of conducting material, that con- 
nects the instrument bo the lunar surface through a voltage 
stepping supply. Figure 2-7 is a schematic representation if 
this connection. The voltage supply is also connected to a 
wire grid above the ion entrance aperturss. This grid is 
at instrument ground . 

The purpose of the ground plane stepper is to compensate 
for any existing Ixinar surface potential and to determine its 
value. For example, if the lunar surface potential is +10 
Volts with 0 volts across the stepping supply and the instru- 
ment is electrically connected to the surface, then a 10 eV 
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Figure 2-4 . Diagram of the internal configuration of 
the various components of the SIDE, [after Freeman et al 
(1969) ] . 

Figure 2-5. Results of mass analyzer calibration, 
(plotted from data compiled by Linderaan [1973]) . 

Figure 2-6. External view of SIDE showing ground plane 
and grid over the entrance aperature of the instrument, 
[from Freeman et al. (1969)]. 

Figure 2-7. Schematic diagram showing relation between 
grid and ground plane. The ground plane is tied to the 
aperature grid through a stepping voltage supply . 
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proton would not be detected due to the +10 V potential at 
the entrance aperture. If, on the other hand, the voltage 
supply maintained a -10 V potential on the grid, with respect 
to the ground plane, the net potential seen by an incoming 
10 eV proton would be zero and it would enter the detector as a 
10 eV particle and would be counted as such. This connection 
between the lunar surface and the instriiment through the step- 
ping voltage supply can also be used to deteonnine the value 
of the lunar surface potential. If we assume that the ambient 
lunar atmosphere consists of essentially thermal ions {energy 
0 eV) then as the grotind plane is stepped through various 
voltages, the ions will begin to appear in the energy channel 
corresponding to the absolute difference between the lunar 
surface potential and the voltage on the grid. For example, 
if the voltage on the grid is — 2C V and particles show up in 
the 10 eV/q channel, then we can expect that if these particles 
were indeed from the ambient lunar atmosphere , having almost 
no energy initially, that the lunar surface potential is about 
+10 V. 

The usefulness of the ground plane stepper in determining 
the lunar surface potential has been demonstrated by Fenner 
et al. [1973J , and it has been determined that a potential of 
about +10 V is appropriate for the dayside of the moon. Based 
on data collected near the terminators a value of about -100 V 
has been set for the night side. Figure 2—8 shows the 
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Figure 2-8 . Results of determination of the lunar 
surface potential in two regions of the orbit. 


IPenner et al. (1973)]. 
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lunar surface potential data thus far obtained. Theoretical 
predictions of this potential, based on various values of 
photoelectric emissivity of the lunar surface, have been made 
by Manka [1972J and these theoretical predictions are also 
displayed in figure 2-8. 


2.4 Fliix measurement 

The counting rate R, which is the quantity actually avail- 
able from SIDE data, can be converted into a more useful quan- 
tity, namely the unidirectional differential energy flux 

j ( ^ , n ) . 

The counting rate at energy ^ is given by, 



g( ^ , 


0 r ) 3 { ^ , 8 , fp ) 6.6 d<p 


2-13 


g C ^ r 8 , 0 ) =unidirectional geometric 
factor at energy ^ 


and 


j { ^ , 0,0 ) =unidirectional differential 
flux 

If we assume that the flux is independent of 8 and <p over 
the solid angle of the detector's field of view so that 
3 { ^ , 8 r ) is approximately constant, then. 


R =y*j ( ^ » 0 , S6) J g{ f , 0 , 9 b)d 0 d0 d^ 


‘'n 

n)G{^)d^ 
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where , 


Gie) 



9i ^ f 6 , 


= H(^) €[^) 
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and 5C^)=A{^) OC^) and e(^)=counter efficiency, 
eff eff 

A is the effective area of the entrance aperture as seen 
eff 

by an incoming particle and is the effective solid angle. 

If we also assume that the flux is constant over the energy 
bandwidth of the detector, then 

R = j G{^)di 

R = j ( f , O ) 2-16 

G = / G ( ^ ) d^ (cm^ ster ev) 

where the integration is carried out over the energy bandwidth. 

Therefore, the unidirectional differential flux at energy 
^ is given by. 


j( ^,n)=R/G^ 


, -2 -1 ^ -1 
(cm sec ster eV ) 
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likewise, the omnidirectional differential energy flux is. 
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Figure 2-2 shows the results of the geometric factor cal- 
ibration for the Apollor XIV instrument [Lindeman, 1973] . 

2.5 Deployment 

The SIDE is part of a group of experiments, known as the 
Apollor Lunar Surface Experiments Package (ALSEP) , designed 
to investigate various phenomena associated with the moon it- 
self and its interaction with the solar wind and the earth s 
magnetosphere. SIDE'S have flown on three Apollo missions; 
all were successfully deployed and all are still operating 

at the time of this writing. 

The Apollo XII site is located at about 3“S latitude 
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and 23 °W longitude in Oceanus Procelarum (Sea of Storms) . 

The Apollo 12 SIDE was placed on the lunar surface by Alan 
Bean and Charles Conrad during the first EVA which began at 
about 6 a.m. EST November 19, 1969. 

The Apollo XIV mission took place in early 1971, result- 
ing in the placement of the Apollo 14 SIDE on the lunar sur- 
face at about 3°S latitude and 17 °W longitude in the Fra Mauro 
formation. Alan Shepard and Ed Mitchell deployed the ALSEP 
on their first EVA which began about 9 a.m. EST, February 5, 
1971. 

The Apollo XV SIDE was deployed by Dave Scott and James 
Irwin in the Hadley-Apennine region of the moon at about 26 °N 
latitude and 3°E longitude again during the first EVA which 
began about 8 a.m. EST July 31, 1971, 

Figure 2-9 shows the positions of the lunar surface of 
the landing sites and their seleonographic coordinates. 
Photographs of the three instruments are shown in figure 2-.1 0 , 
The Apollo 12 and 14 SIDE'S were equipped with a bubble 
leveling device attached to the top of the instrument that 
allowed the astronauts to align the vertical axis of the SIDE 
with the local vertical to within a 5“ accuracy so that the 
look direction of the detector is approximately in the eclip- 
tic plane. The Apollo 15 instrument was equipped with an ad- 
justable leg that allowed compensation for the high latitude 


T 
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Fi gure 2-9. Seleonographic locations of Apollo 12, 

14 and 15 ALSEP sites. 

Figure 2-10 . From left clockwise, photographs of 
Apollo 12 and 14 SIDE'S deployed on the lunar surface, 
Apollo 15 SIDE in the lab. [after Medrano, 1973]. 
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deployment so that its detector look direction is also approx- 
imately in the plane of the ecliptic, (see figure 2-lla) 

The look direction of the actual detector in each of the 
instruments is at a 15° angle from the instrument x— axis, which 
coincides with the local vertical at the ALSEP site for the 12 
and 14 instruments. This 15° offset determines the actual look 
direction, in longitude, with respect to the moon. The config- 
uration in longitude is illustrated in figure 2-llb. Figure 
2-12 shows the body axes of the three instruments. The relat- 
ionship of the body axes to the local coordinate system at 
ALSEP is as follows. 

The local coordinate system is known as system 1, and in 
the discussion in chapter 4, it is referred to as S' ' ' . It is 
defined by an x-axis that is aligned with the local vertical 
at ALSEP, a z-axis that is perpendicular to the x-axis and dir- 
ected northward along the meridian of deployment. The y— axis 
completes the orthogonal right handed system. That is, y 
points east. 

In the case of the Apollo 12 and 14 SIDE'S the body x-ax- 
is of the instrument coincides with the x-axis of system 1. 

The body x-axis of the 15 instrument is in the plane of the 
meridian of deployment but tilted southward by 26°. The body 
y-axis of the 15 instrument coincides with the y-axis of sys- 
tem 1 and the z-axis is orthogonal to x and y. The body z-ax- 
is of the 14 instrument coincides with system 1, but the body 
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Figure 2-11. 


Figure 2-11.. 


Figure 2-12. 


a. shows the configuration in latitude 
of Apollo's 14 and 15. Apollo 12 is ap- 
proximately represented by the drawing 
for Apollo 14 . 

shows the configuration in longitude 
of all three SIDE’S. Apollo 14 SIDE looks 
most nearly tov/ard earth. 

F.elationship of the instrument body axes 
to the local ALSEP coordinate system 1. 
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y- and z-axes of the 12 SIDE are the negatives of system 1. 

For further discussion of the relation of the body axes 
to other relevant coordinate systems, see chapter 4 and Hills 
[1969] . 

The deployment of the instruments in this fashion allows 
a directional sampling of magnetospheric and solar wind part- 
icles at any given time. The look directions of the three de- 
tectors for various positions in the lunar orbit are shown in 
figure 2-13. Note that the Apollo 14 SIDE looks almost direct- 
ly back at the earth at any position in the orbit, which makes 
this instrument particularly suitable for the observation of 
particles streaming outward from the earth's bow shock. Be- 
cause of the excessively high background counting rate in the 
Apollo 15 instrument and the discontinuous nature of the Apollo 
12 data for the period of interest in this thesis, only data 
from the Apollo 14 SIDE were used. 

2.6 Magnetometer [Colburn, 1969] 

The data on the interplanetary magnetic field and those 
from the geomagnetic field were collected by the NASA Ames mag- 
netometer onboard the Explorer 35 satellite which achieved lun- 
ar orbit on July 17, 1967 with a periselene of 1.4 (lunar 
radii) , an aposelene of 5.4 and an orbital period of 11.5 
hours . 

The magnetometer contains three orthogonal fluxgate sen- 


Figure 2-13. Longitudinal look directions of the 
three SIDE'S for various positions in lunar orbit. 
Note that the Apollo 14 SIDE looks most nearly toward 
earth. Different sunset and sunrise times are due to 
the difference in longitude of deployment. 


[Hills et al . (1971) ] . 
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sors which measure the ambient magnetic field in the neighbor- 
hood of the spacecraft and a vector measurement of the field 
is obtained every 6.14 seconds in the 20 gamma range. The 
signals from the sensors are spin-demodulated onboard to pro- 
vide a vector in a fixed spacecraft coordinate system, V7ith 
the final data available from Ames in Solar Equatorial or Solar 
Magnetospheric coordinates, averaged over 81.8 second intervals 
The data are further processed on the SDS 910 computer at 
Rice to provide the field components in Solar Ecliptic coordi- 
nates , [Lindeman ,1972] . 

The solar ecliptic coordinate system is defined as follows 
The x-axis is defined by the earth-sun line with positive x 
toward the sun. The z-axis is normal to the plane of the ec- 
liptic and positive z is toward north. North is the direction 
defined by the direction of the orbital angular momentem vec- 
tor of the earth, y completes the right-handed system. Solar 
ecliptic longitude is defined as 

tan“^ (y/xj 

and solar ecliptic latitude is defined as 

Q= tan~^ (z/(x^ + y^) 2-22 

Due to a failure in the sun sensor in the spinning space- 
craft, the values of the solar ecliptic longitude are often 
in error. It is possible, however, to get a rough correction 
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for this error by observing the longitude data in the geomag- 
netic tail region where stable neutral sheet crossings are ob- 
served. 

If a neutral sheet crossing occurs when the spacecraft is 
near 180° in longitude, one would expect that while in the 
southern lobe of the field, the solar ecliptic longitude of 
the field would be about 180° and after the crossing, the 
field would have a longitude of approxiamtely 0°. By observ- 
ing the actual measured values and know^ing the nominal values 
of the field longitude, the correction angle can be determined. 
An example of this is shown in figure 2-14 where a correction 
angle of about 28° was obtained for the lunation with which 
this thesis is concerned. For the present analyses, the un- 
corrected values of the longitude were used. Any physical 
implications of using these uncorrected longitudes are dis- 
cussed in chapter 5 . It should also be noted that the Explor- 
er 35 magnetometer has no onboard recorder and therefore there 
are frequent data gaps of about 5 hours during each orbit due 
to the passage of the satellite behind the moon. 

2 . 7 Other lunar surface experiments 

The SIDE is only one of a variety of instruments design- 
ed to investigate plasma properties in the vicinity of the 
moon. These other lunar surface experiments provide data 
which can potentially be correlated with SIDE observations . 
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Figure 2-14 ■ Solar ecliptic longitude of the magnetic 
field near the anti-solar position in lunar orbit. The 
technique for determining the correction angle is illus- 


trated. 
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Although detailed discussions of chese other experiments are 
beyond the scope of this thesis, the descriptions below are 
included to provide an insight into the flavor of the particles 
and fields research being conducted on the lunar surface. 

2 . 8 CPTpT^E-Charged Particle Lunar Environment Experiment 
The CPLEE consists of two charged particle analyzers, 

each containing five C-shaped channeltrons as negatively 
charged particle detectors and one channeltron as a positively 
charged particle detector, with electrostatic deflection plates 
providing the energy discrimination capability. Particles 
with energies from 50 eV to 50 keV can be detected by CPLEE 
from two directions simultaneously; one analyzer looks verti- 
cally with a 4°X20° field of view and the other detector is 
directed at 60° from the local vertical with the same field 
of view. 

The CPLEE has observed such phenomena as a photoelectron 
layer at the lunar surface, solar wind particles and penetra- 
tion of particles of auroral energies to lunar distances in 
the geomagnetic tail [O'Brien, 1971]. 

2.9 CCGE-Cold Cathode Gauge Experiment [Johnson et al., 1971] 
The cold cathode gauge experiment is designed to measure 

the neutral gas pressure of the lunar atmosphere. The instru- 
ment is attached to the SIDE and is connected to the ALSEP 
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central station through the SIDE, but it is a separate exper- 
iment consisting of a chamber in which a spool shaped cathode 
is connected to an electrometer. The cathode is surrounded 
by a cylindrical anode that is held at +4500 volts. When 
neutral gas enters the chamber, a self sustaining discharge 
occurs with the ion current through the cathode being a measure 
of the gas pressure. 

Results from the CCGE set neutral gas number densities 

5-3 7-3 

for lunar nightime of about 10 cm and about 10 cm for 
daytime with the daytime value partly due to local experiment 
outgassing. 

2.10 LSM--LPM-Lunar Surface Magnetometer and Lunar Portable 
Magnetometer . iDyal et al., 1972] 

The liinar surface magnetometers each consist of three 
orthogonal fluxgate sensors, mounted on booms, with the elec- 
tronics package located at the junction of the booms. This 
configuration provides a vector measurement of the static 
as v/ell as the time varying magnetic fields at the ALSEP sites . 

The lunar portable magnetometers are self contained and 
also provide vector measurements of the field at several dif- 
ferent locations along the astronaut traverses. Table 2-2 
shov 7 S the range and sensitivity of the Apollo 16 LSM and LPM 
[Dyal et al., 1972]. Several of these magnetometers have 
flown on Apollo missions and a summary of the preliminary 
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Suiitimry of Lumr Surface Magnciic Field Measitremenls 



Site 

CoordiMalM, deg 

Field magnitude, y 

Magnetic field coniponerirj, 7 

^up 

^cast 

^north 

Apollo If’ ^ 
Al.SliPsite 

B.9 S 15.5 E 

231 

-174 

-73 

134 

Site 2 


180 

-179 

-10 

19 

Site S 
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-25 

-25 

Situ 13 
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-141 

-209 

-186 
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-42 

-93 

66 

Apollo 15 A1.SEP site 

26.1 N 3.7 E 

6 +4 

4+4 ! 

1 +3 

4 ±3 

Apollo 14 land big site 
Site A 

3.7 S 17.5 W 

103+5 

-93 +4 

38 +5 

-24 +5 

Site (' 


43 + 6 

-15 ±4 

-36 +5 

-I9±8 

Apollo 12 Al.SliP site 

3.2 S 23,4 W 

38+2 

-25.8 + 1.0 

11.9 +0.9 

-25.8 ± 0.4 


“•riiL> error in ;iU component:; ol Die Apollo U> rielil r.-:ltliiu;s is +25 7 . These limits will be redticctl al'lcr Imther data 
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Estimated ^Ue i-lux A vcrcii;cs During the Times 


of Foil Exposure 


Mission 

Exposure inUiation 

Exposure dura/ton, 
/U’smitt 

A verage soiar wind 
‘^ite flux. 
l(F‘cin'~sce'^ 

Dale 

O.III.I., 

Iir.'mi/i 

Apollo 11 

July 21. 1969 

03:35 

01:17 

6.2 1 1.2 

Apollo 12 

Nov. 19, 1969 

12:35 

18:42 

8.1 J 1.0 

Apollo 14 

Ixb. ,S, 1971 

15:15 

21:00 

4.2 i 0.8 

Apollo IS 

July 31, 1971 

19:36 

41:08 

17.7 1 2.5 ■ 

Apollo 16 

Apr. 21, 1972 

23:01 

45:05 

12.0 ! 1.8 


Solar Wind Abundance Ratios Determined From 
the Foil Exposure Periods of the Apollo Missions 


Mission 


‘>Hcl~^Ne 

-f>Ne/--Ne 

Ne 

Apollo 11 
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- 

Apollo 15 
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13.80 ±0.40 
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Preliminary Results on Argon 
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First Results From the Analyses of the A luminum Foil 
From the Apollo 16 5ilT Experiment 
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results of ambient field measurements at different locations 
is shown in table 2-3 . 

2.11 SWCE-Solar Wind Composition Experiment iGeiss et al., 

1972] 

The SWCE is one of the simpler lunar surface experiments 
and consists of exposing an aluminum foil strip, 30X140 cm, 
directly to the solar wind and then returning the foil to 
the laboratory for analysis . The foil is baked out in a 
vacuum chamber and in some cases melted to release the trapped 
Solar wind particles which are subsequently analyzed with a 
mass spectrometer. The SWCE was one of the experiments flown 
on the Apollo 11 mission as well as several others and a sum- 
mary of the results from these missions is shown in Table 2-4. 

2.12 SWS-Solar Wind Spectrometer iSynder et al,, 1970] 

The solar wind spectrometer was designed to investigate 
properties of the solar wind and its interaction with the moon. 
It consists of seven symmetrically arranged Faraday cups, one 
directed vertically, and the other six at 60° angles from the 
vertical. This configuration allows observation of any solar 
wind anisotropies that may exist. Energy spectra may be 
obtained by modulating the voltage on the grids that precede 
the Faraday cups and noting the cup current changes. The 
SWS cycles through 14 energy steps ranging from 10 eV to 
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10.4 keV for protons and from 10 eV to 1480 eV for electrons 
in the high energy range and each energy step is reduced by 
a factor of 1.68 in the low energy range. 

Solar wind proton spectra obtained in this way show pro- 
ton velocities of typically 400 to 600 km/sec. The SWS does 
not generally see solar wind particles during the lunar 
night, indicating a plasma void behind the moon. 

The power for all of the experiments described, with the 
exception of the LPM and the SWCE is provided by a Radioisotope 
Thermoelectric Generator (RTG) which provides a continuous 
70 W. The individual experiments are connected to the ALSEP 
Central Station by Mylar coated ribbon-like cables where the 
data from each of them are collected and transmitted to earth 
through a high gain antenna mounted atop the central station 
and aimed at the earth. These data are received at various 
earth based tracking stations, recorded on magnetic tape and 
distributed by NASA to the various investigators. 
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Chapter 3 

This chapter is a description of the model used to study 
the behavior of suprathermal protons that are ejected into 
the solar wind in an upstream direction fron> the earth's bow 

shock . 

3.1 Equations of motion 

Consider a particle being ejected from the earth's bow 
shock with some initial velocity and define a coordinate 
system in which the z-direction is parallel to the interplan- 
etary magnetic field B, y is in the direction of ^ ® 

and X completes the right-handed system. If the origin of 
this coordinate system is fixed at som.e point on the surface 
of the shock, then the particle sees an electric field 

^ = -V X"B, due to the convection of the interplanetary 
sw _i 

field lines past the earth at the solar wind velocity 

The equation of motion that describes the particle s 

trajectory upon leaving the shock is, 

P = q( E + V X B ) 3-1 

dt P 

In a coordinate system V 7 ith B = Bz, one can choose E = Ey 
without loss of generality since E is always perpendicular to 
"b. The solution of equation 3-1 gives the time dependent 


50 


coordinates of the particle trajectory. 


x(t) = 


V 


V -V „ \ V 

^^ 7 -; sinCOt + ( 1 -cosCiJt) + v t 

UJ I CO d 


3~2 


V 

y (t) = sinCdt + 

03 


V -V, 

xo d 
03 


(cosOJt-1) 


3-3 


z (t) = V t 

zo 


3-4 


Equations 3-2 and 3-3 are the parametric equations of a 
cycloid in the x-y plane. For a derivation and discussion 
of these equations , see the appendices. The motion in the 
z-direction is uncoupled to motion in the x-y plane and hence 
if there is no initial velocity in the z-direction, the part- 
icle will move entirely in the x-y plane. 

If we transform the above equations to a frame of reference 
that is moving parallel to the x-direction with velocity 
we see that the particle describes circular motion with center, 


mv 


c qB 


yo 


and radius 


R = (m/qB) 


Y =-2^(v--v ) 

c qB d xo 


2 2 
(v -v.) + V 

xo d yo 


3-5 


1/2 


3-6 


in the moving reference frame. That is, the guiding center 
of the circle moves in the x-direction with velocity V^=E/B, 
the drift velocity, (see figure 3-1) 
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Figure 3-la. Illustration of particle drift motion 
for a proton' in crossed E and B fields when the initial 
velocity of the particle is zero. 

Figure 3-lb. View of drift motion in the plane 

showing the effect of adding a finite z-component of 
the velocity. 

Figure 3-lc. View in the plane showing the 

particle's resultant guiding center trajectory. 
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Due to the uncoupling of the z -direction motion from the 
x-y planar motion, a simple expression can be v/ritten that 
determines the resultant motion of the guiding center of a 
P^^rticle that leaves the earth's bow shock with 
some intial velocity V^. The resultant velocity of the 
particle's guiding center is given by. 


V 

r 




3-7 


where the component of the particle's initial velocity 

that is parallel to the interplanetary field and is given 
by. 



= (E X B) /B^ 


3-8 


3 • 2 Guiding center behavior 

Figure 3-2 is an idealized representation of the magnet- 
ospheric cavity as it exists in the free flowing solar wind. 

In a solar ecliptic coordinate system the interplanetary mag- 
netic field has a solar ecliptic longitude of approximately 
315°, the "garden hose angle", and lies more of less in the 
ecliptic plane. 

To understand the basic longitudinal behavior of the model, 


Figure 3-2. Diagram demonstrating the asymmetry of 
the ^ X B drift model and the basic guiding center 
behavior of a zero pitch angle particle. 
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consider a particle of zero pitch angle leaving the shock at 
point P in figure 3-2 and traveling along the field lino with 
a velocity V ^ determined by its initial energy. If the in- 
terplanetary field were stationary, the particle would move 
along the path determined by and it would intersect the 

lunar orbit at point B. Due to the convection of the field 

* 

lines past the magnetospheric cavity at the solar wind veloc- 
ty, there exists an interplanetary electric field E= B 
such that the particle's guiding center follows a path that is 
detearmined by, 


V = V + V 
r sw 



B 


3-9 


and the particle will actually follcw a guiding center path 
that intersects the lunar orbit at point A. If the moon is 
also at point A and the detector is looking in the right dir- 
ection, then the particle can be seen there if it is of the 
proper energy to be detected by the SIDE . 

Clearly, if a particle leaves point P with less energy or 
a higher pitch angle, i.e, a "shorter” V , the same drift vel- 

ocity vector added to its V will yield a resultant velocity 

7 / 

vector that defines a guiding center path that intersects the 
lunar orbit at a point that is closer to the shock than point 
A. Conversely, a particle leaving the shock with greater par- 
allel energy, i.e. a "longer" 


, v;ill have a resultant 
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velocity vector that will define a guiding center path that 
intersects the lunar orbit at a point further away from the 
shock than point A. 

That is, whether or not a particle of a given energy 
is detected at the lunar orbit is dependent upon the orienta- 
tion of the interplanetary B field, the magnitude and direc- 
tion of the solar wind and the particle's initial energy and 
pitch angle upon leaving the shock. 

From the above description of the energy selection mech- 
anism that is operating in the post shock region of the 
lunar orbit, it becomes clear that for a particle leaving 
the shock with a given parallel energy E. , and hence a given 




initial parallel velocity Vp , that there exists a point Q 


on 


ii 


the bow shock that defines an energy cutoff point in the lunar 

orbit. If the moon has advanced in its orbit beyond the point 

C in figure 3-2, then particles of energy E-n leaving the 

nl 

bow shock can no longei be seen by the detector. 

The cutoff point varies with energy and interplanetary 
field conditions, but one expects that in the nost shock cros- 


sing region of the lunar orbit, that the lower energy particles 
will begin to disappear from the data as the moon advances 
toward instrument nightime. This effect is observed in the 
data and will be discussed again in the chapter concerning 
the data, (chapter 5) 
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A completely analogous mechanism is operating on the inbound 
side of the bow shock and the appropriate vector diagram for 
that region is also illustrated in figure 3-2. Due to the 
geometric asymmetry in the field configuration, the data from 
the inbound side of the shock will show a somewhat different 
effect. The higher energy particles will be seen further away 
from the shock than v/ill the lower energy ones, but none will 
be seen as far away from the shock as on the outbound side for 
a given field configuration. This is in agreement with the 
"dusk-dawn" asymmetry observed by Asbridge et al. [1968]. 

The above discussion required that the solar ecliptic lat- 
itude of the field be approximately zero since the moon's 
orbit lies approximately, (+ or - 5°) , in the ecliptic plane. 
This would suggest that if the interplanetary B field were 
not in the ecliptic plane or close to it, then particles leav- 
ing the bow shock might not be seen by the detector. Indeed, 
if the field is far away from the ecliptic plane, particles 
leaving the shock along these lines could likely miss the 
lunar orbit entirely, (see figure 3-3) 

In the above discussion, the interplanetary field was 
taken to be approximately in the ecliptic plane. If the 
solar wind vector also lies approximately in the ecliptic 
plane, this requires that E=- B be perpendicular to the 

ecliptic plane and hence will also lie in the ecliptic 
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Figure 3-3. Illustration in three planes of the dep 
endence of the resultant guiding center trajectory on 
the interplanetary field latitude. 




I 
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Figure 3-3 
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plane as will the resultant velocity vector % . But if B 
is not in the ecliptic, the resultant velocity vector will 
be out of the ecliptic plane. As can be seen in figure 3-3, 
whether or not a particles 's guiding center trajectory inter- 
sects the lunar orbit is fairly sensitive to the interplanet- 
ary field latitude. An analysis of the data shows this effect 

and will he discussed in chapter 5. 

It should also be pointed out that there is an interesting 

dusk-dawn asymmetry in the direction of the interplanetary 
"E field when the 't field is at high positive or high negative 
latitudes. Figure 3-4 shows that for a B field latitude of 
+90° on the dawn side, that the interplanetary electric field 
points outwardly from the shock and similarly a B field lat- 
itude of -90° on the dusk side generates an E field that also 
points away from the shock. This gives rise to another mech- 
anism, as discussed by Manka [1972], for getting particles 
into the detector that has a definite dusk-dawn asymmetry 
associated with it. 

The foregoing description of the model has ignored any 
influence of the perpendicular energy of the particle and 
hence the size of its gyroradius . To gain understanding of 
this influence, the actual particle trajectories will now 

be discussed. 
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Figure 3-4. Diagram illustrating the dusk-dawn asym 
metric interplanetary electric field E that occurs at 
high positive and negative values of the interplanetary 


B field latitude. 
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3 . 3 Actual particle trajectories 

The gyroradius of a 3 keV proton in a ST field is about 
1500 km which is comparable to the lunar radius (-^JlTOO km). 
This implies that a strict interpretation of the above treat- 
ment, where the particle follows the guiding center path det- 
ermined by is not quite correct. The actual trajectories 

depend on the particle’s initial velocity and pitch angle a 
and can be calculated from the parametric equations 3-2, 3-3, 
and 3-4. The pitch angle is the angle between the initial 
velocity vector and the direction of the field which implies 
that 



V = V cos a 
20 o 


3-10 


where 


and 


O \ XO yo 20/ 


(2 (f/m) 


1/2 3-11 


V - V sxn a = 

1 . 


4- v2 ^ 
XO yo 
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For typical conditions, Vg^p400 km/sec, B=57 , several 
cases of proton motion are shown, to scale, in figures 
3-5 through 3-9. It is evident from these figures that part- 
icles whose guiding center trajectories intersect the moon 
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Figure 3-5. Scaled illustration of cycloidal motion 
of typical bow shock protons. Top curve is the case 
of all parallel energy and the bottom curves are cases 
of 1 keV of perpendicular energy distributed in three 
different ways. and are the values of the part- 

icle's initial velocity in the x and y directions res- 
pectively in units of km/sec. 

Fiqu3^6S 3-6 through 3-9 . Additional cases of cycloidal 
motion for perpendicular energies higher than 1 keV. 

Note the interesting case in figure 3-6 where the x-comp- 
onent of the initial velocity is approximately equal to 
the drift velocity (283 km/sec) . (See equation 3-6) 
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can enter the detector from a wide range of angles. The 
implications of this on the interpretation of the data are 
discussed in chapter 5 . 

3.4 Asbridqe Model 

The model proposed by Asbridge et al. 11968] is very 
similar to the present treatment and the basic idea behind 
their model is illustrated in figure 1-6 . The resultant 
velocity vector of the particle is determined by / 

where/Vj_ is the particle’s initial velocity parallel to 
the field. This differs from the expression for the re- 
sultant velocity vector given in equation 3-7 by the factor, 

- { V • B ) ^ 3-14 

sw 

which is just the component of the solar wind velocity par- 
allel to the field. WhenV^^^is' perpendicular to B; the two 
expressions are equivalent . 

The basic difficulty with eqation 3-13 is that it appears 
to give the resultant guiding center velqcitsr veqtpr^^^^;^^ 
frame of the bow shock for a particle that is injected into 
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the solar wind along the B field in the shock frame. There 
is no question that if a particle left the bow shock in the 
direction defined by that it would continue along this 
path undeviated by any electric fields since 

■^ = ( ^ + Vg X B ) = 0 

nonrelativistically . However, if the particle actually had 
a velocity in the shock frame as indicated in figure 1“6 , 
it would not go in the direction defined by as indicated 

by the Asbridge model. 

point can be clarifxed by considering the transform 
atioh between a frame of reference that is moving with tbe 
solar wind and the rest frame of the bow shock. Since Vg^«c, 
the transformation is Galilean, 

X = X' + Vg^t 

y = y' 3-16 

2 = 

where the moving frame is being translated along the x-axis at 
the solar wind velocity. This implies that. 
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"v ="v’ + V • 3-17 

sw 

If V' is directed along "B in the moving frame, then the 
velocity seen in the rest frame is V. Equation 3— 13 is actual 
ly just a transformation from a frame moving with the solar 
wind to the rest frame of the shock, i.e., 

3-18 
3-19 

It is not realistic to assume that one knows the inital 
velocity of the particle in the moving frame and therefore 
equation 3-7, derived in this chapter, is the appropriate 
expression to determine the resultant guiding center trajec- 
tory of a particle in the rest frame of the shock. It is 

-X 

interesting to note that when is not perpendicular to B, 
equation 3-9 yields a resultant guiding center trajectory 
that allows particles to travel to points at the lunar orbit 
that are more sunward than those predicted by equation 3-13 , 

3,5 Conclusions 

The discussions in this chapter lead to the following 
conclusions, 


and 


^i 


= V 


= V 
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a. ) There are particlular interplanetary field longitudes 
that are favorable to the observation of particles of a given 
energy at a given point in the lunar orbit and this fact gives 
rise to energy cut-offs in the post shock region of the lunar 
orbit and energy cut-ons in the pre shock region. 

b. ) There are particular interplanetary field latitudes 
' that are favorable to the observation of particles in the 

post and pre shock regions of the lunar orbit, viz. those in 
the neighborhood of zero degrees. 

c. ) There is an inherent dusk-dawn asyininetry in the model 
that should be observable in the data. 

d. ) There is a dusk-dawn asymmetric interplanetary elec- 
tric field mechanism for getting particles into the detector 
that is dependent on high positive and high negative values 
of the interplanetary B field latitude. 

e. ) Although the Asbridge model of guiding center tra- 
jectories gives the correct qualitative picture of particles 

<» 

leaving the bow shock and the correct quantitative picture if 
interpreted correctly, the E X B model proposed in this thesis 
is more appropriate for the reasons discussed in this chapter. 
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Chapter 4 

-In light of the foregoing discussion, a computer model 
for the "e X "b mechanism has been generated. The description 
of this computer model and its connection to the physical 
model of the previous chapter will be given below. In order 
to generate a qualitative picture of V7hat the detector should 
see as a function of time, it is necessa3:y to calculate the 
components of a vector '^(t) that represents the look direction 
of the detector as a function of time and use this vector in 
conjunction with the vectors V^( ^ , Q! ) , the resultant guiding 
center trajectory vector/ and r/ the position vector of the 
moon to determine whether or not a particle is seen. By de- 
fining a source with certain characteristics and "acceptance 
criteria" for the detector, a spectrum of counts vs. energy 
can be obtained and compared to the real data. What follows 
is a description of the computer model. 

4.1 Orbit 

In order to determine whether or not a particle leaving 
the bow shock will hit the moon, it is necessary to determine 
the position vector of the mpon, /r. Figure 4-1 shows the or- 
ientation of the lunar orbit with respect to an inertial sys- 
tem whose x—y plane is in the plans of the ecliptic. The or- 
bital inclination , B , is exaggerated in the figure for clarity 
it has a value of about 5°. The x’ axis intersects the lunar 
orbit at perigee (denoted by tt ) in the figure . 
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Figure 4~1. Illustration of the relation of the lunar 
orbit to the solar ecliptic plane showing time depend- 
ent angles relevant to the calculation of the position 
of the moon . rt denotes perigee . Q is exaggerated in 
the figure for clarity. In reality Q A-* 5°. 



Figure 4-1 
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and precesses with a period of approximately 8 . 5 years in a 

counterclockwise direction when viewed from a point on the 

z> -axis above the x^-y' plane. The line that divides the time 

between when the moon is above the ecliptic plane and the 

time when it is belox '7 the ecliptic plane is knoTO as the line 

> 

of nodes and it precesses in a clockwise direction when look- 
ing in the negative z-direction from a point above the eclip- 
tic plane. The line of nodes precesses with a period of ap- 
proximately 18.6 years. 

The moon moves in a slightly elliptical ( €>^.05) orbit 
that lies in the x'-y’ plane with perigee occuring on the 
x’ axis. The distance from the earth to the moon as a func- 
tion of time is given by , 


a (1- 6^ 

1 + e cos ( w t- t|j ) 


where, a = semi-major axis, €= eccentricity, W= angular 
velocity of the moon based on the sidereal period, argu- 

ment of perigee, and t begins at ascending node. 

The angles, relevant to the Eulerian transformation, as 
functions of time are determined by , 

,^(t) = ~ %on^ implies CW precession 


(p ~ argument of line of nodes at t 0 


^lon" ^ "^lon 


T = 18^.6 yrs, 
Ion 
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t (+) implies CCW precession 

« loa 


= argument of perigee at t-0 


4-3 


ft, = 27T/ T. 
loa loa 


T. =8.6 yrs. 
loa 


0 (t) = const. 




The transformation from a coordinate system, S*', that is 
fixed in the moon , to the solar ecliptic system , is accom- 

plished by first transforming from S'* to S ’ and then from 
S' to S and finally from S to S^^. That is if ^ represents 
the matrix transformation from S'* to S ' , ^represents the 
transfo 3 rmation from S ' to S and ,S_ represents the transforma- 
tion from S to S^^ then, 

"r ' = ,A. ~r' ' 

r =_B r' 4-4 

r = ^ r 
se — 

which implies that, 

=r S B A ■? ' * 
se 

The coordinate system S'* is simply rotatxng with respect 
to S* so that the angle between the x* '-axis and the x*-axis 
is, 

0= &j t- t/f(t) 4-5 
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and since z ' ' ?= z ' , 

A = 


/ 


cos 0 


-sin© 



sin© COS0 0 



0 
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The transformation from S' to S is given by equation 4-7, 
shown on the next page. 

The inertial system S is stationary with .respect to the 


fixed stars but the x-y axes of the SE system rotate slowly 

with respect to the inertial system due to motion of the earth 

in its orbit about the sun. That is, the -axis rotates with 
respect to the x-axis about a degree per day in the counter- 
clockwise direction as viev/ed from above the ecliptic plane. 
Therefore, the transformation from the inertial system to the 
solar ecliptic system is given by the matrix, 

S = 

where , 

i= io * ’■e 

T = sidereal period of the earth 

e 

Therefore the components of the vector r in the SE system 


/ 


cos ^ 
sin ^ 




-sin i 
cos i 
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are^ 


where , 





C S B A 
( S B A 
( S B A 


1 t 


'X 



1 I 
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and the matrix elements of SBA are shown on the next page. 


4.2 Detector Look Direction 

If one know the orientation of the detector look vector, 
"d* * ' in a coordinate system, S''", that is fixed at the ALSEP 
site f then the components of this vector can be expressed in 
a S.E. system, S , by making an Euler ian transformation 
from S’’’ to Sgg. 

The moon's rotational axis makes an angle with the 
ecliptic plane of about 1.5 degrees and the rotational period 
is equal to the sidereal period. Figure 4-2 shows the rela- 
tionship between the coordinate system fixed at the ALSEP 
site and the solar ecliptic system. , 

The Apollo 12 and 14 instruments were deployed very near 
the lunar equator and for the purposes of this treatment we 


cos f L ® (cos y}f cos<^ - cos Q sin ^ sin ^ ) • 
-sin j^cos 0 (cos^ sin0 + cos 0 cos^ sin 0^) ■ 

sin^ l^cos © (cos 0 COS0 - cos0 sin0sin0)- 
+COS ^ ■ jcos 0 (cos0 sin0 + cos Q cos 0 sin 0 ) 

(cos© sin0 + sin© cos 0 ) sin0 

Equations 4-9 




sin© (sin0 cos0 + cos0sin0cos0)J 
sin 0 (sin0 sin0 - cos 0 cos0 cos 0 ) 

sin© (sin0cos0 + cos0 sin0cos0) 
•sin© (sin0 sin0 -cos0 cos0 cos 0)1 


i 




u> 


Figure 4-2. Relationship of the ALSEP coordinate sys- 
tern S*'', or system 1, to the SE coordinate system. For 
deployment at the sub-earth point, x'*' points to earth 
and remains approximately that way due to the equality 
of the lunar sidereal period and the rotational period,. 
Small changes are due to physical librations . 



const 

X«)=X.+3i< 

%4 = sidereal period of moon 
^tc= sidereal period of earth 
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will as Slime that they Were flep iQy ed exact ly on. the egiiator , 
Then at the ALSEp site the unit vector x’ ' ’ is in the direc— 
tion of the local vertical, ^ points north and in the case 
of eguatorial deployment, z’’^ is parallel to the rotational 
axis of the moon, y’ ' ’ completes the right-handed system 
and is directed eastward, SM* - system 1, 

C assume Apollo 14) , was deployed at the 
ALSEP .site so that the body axes were coincident with x'' ■ , 

■y’ ’ ’ and , The vector that represents the look direction 
of the detector in the S ’ system is then, 

i 

D cos 15 
*D' ' ’ - D sin 15 

“ 

The Euler ian transformation from S''' to S is given on the 

following page. Using this transformation D in the SE 

system is , 
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At this point> it is possible to calculate the result 
tant velocity of particles leaving the bow shock, the pos'- 
itioh vector Of the room r, and the detector look direction 
all in a solar ecliptic coordinate systeni, 

4.3 Shape of the shock’ surface ; 

The source of the particles in this computer model is an 
idealized mathematical representation of the magnetospheric 
cavity that emits particles with a flat energy spectrum, 
isotropically at each point on its surface. The mathematical 
shape is a paraboloid of revolution that is fitted to the 
trne shape of the bow shock surface. This fitting procedure 
will be discussed below. The model retains all of the 
qualitative features of the real bow shock and at least two 
of the quantitative features, the distance of the nose to the 
center of the earth and the distance from the center of the 
earth to the intersection of the lunar orbit with the shock. 
The paraboloid is shown in figure 4~3 and the eguation in 
solar ecliptic coordinates is, 


X = + z^ ) + d 4-14 



where d is the distance from the nose of the shock to the cen- 


ter of the earth and D is the distance from the center of the 


Figure 4-3^ Idealized bow shock surface represented 

2 2 2 2 

by the equation X =-d/D (Y +Z )+d in a SE coordinate 
system. 


1 
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Figure 
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earth to the Shock, measured along the y-axis . The shock 
stand-off distance d, can be calculated by equating the Solar 
wind kinetic ^pressure to the magnetic pressure provided by 
the earth’s dipole field so that 
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and the parameter D, can be determined from knowledge of the 
positions of actual bow shock crossings, cf. fig. 4-4a. 
That is by knowing the time of the bow shock crossing, the 
S*E. longitude can he obtained from the ephemeris data. The 
time of crossing can be determined accurately from the SIDE 
data. An example of this is shown in figure 4 -4b. Solving 
equation 4-14 for D and neglecting any z dependence. 



- 2 . 2 
-d r sin 

r cos - d 
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D v/as calculated in this manner, but a constant value for d, 

f/s/12 R ) was chosen over using the expression 4-15. 
e 



Figure 4-4a. Diagram illustrating the technique for 
fitting the idealized shock surface to the observed one. 
Knowing the distance d and measuring and allows 
determination of D from equation 4 ”14. 

Figure 4"4h. Illustration of hp^ the postion of the 
bow shock is accurately determined for use in calculat- 
ing D as discussed above . 
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4,4 ■ Particle ehix^ss'i'on 

Each point of the surface of the shock is taken to be an 
isotropic radiator of particles with a flat energy spectrum. 
That iSr the same number of particles are emitted for each 
energy and the spectrum looks the same in all directions from 
each point of the surface. The total number of particles 
emitted at a given pitch angle a is proportional to the area 
on e unit sphere between"”^^ and-|-^^ , See figure 

The number of particles available for detection is there- 
fore weighted toward higher pitch angles by a factor of sin a . 
If n is the nmnber of particles per square centimeter passing 
through a unit sphere centered at the point/ then the number, 
dK, emitted at a given pitch angle a is, 

dN = ndA 

dN = 2tt n sine da 
dN = const, sin a 

for a constant interval da , 


The total number N , emitted over the hemisphere, center 
ed at each point, is, 



Figure 4-5a. Shows that for an isotropic radiator 
the number of particles emitted at pitch angle o< 
depends on • . 

Figure 4 -5b. Illustrates that all of the source points 
in the field of view of the detector are tested con- 
stituing a nmnerical integration over the source to 
determine the counts observed at the detector . 



Figure 4-5 
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In tlie computer model , n was chosen to be of the same order 
of magnitude as observed in the data and pitch angle inter- 
vals, Aa f of 10° were used. The initial values of the ener- 
gies ranged from 10 eV to 3500 eV. 

4.5 Acceptance Criteria 

There are essentially two criteria to consider in deter- 
mining whether or not a particle is to be counted by the de- 
tector on the lunar surface. 

1 . ) the particle must come from the direction in which 
the detector is looking and 

2 . ) the particle must have a source point on the sur- 
face of the shock, 

* (it is possible to relax this criterion under certain cir- 
cumstances; see discussion in chapter 5i. Due to the finite, 
size of the detector look cone, there is more than one dir- 
ection from which the particles can come and still enter the 
detector. If S is the angle associated with the look cone, 
then sxnce determines the direction of the guiding center 
path and determines the look direction of the detector, 
a particle can enter the detector if the following relation 
is satisfied. 

A A 

D COS § 


4-19 
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The vector ^ detonuines the direction in space of the 

guiding center trajectory and the vector “r determines the 

point in space through which the par’i^icle must pass. If the 

line, determined by V and the point "r, intersects the sur- 

face of the shock, then a source point exists. The equation 

of a line in three dimensions, passing through the point r, 

with direction V , is, 

r 


x-r y~r z-r 
— X = -i- y = - — — -z 
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V 


. V 


V 


rx ry r z 
The simultaneous solution of this equation and the equation 

of the paraboloid of revolution eq. 4-14, yields a quadratic 
in any of the variables x,y, or z. Then if the discriminant, 
C " 4 A C is real, that is if 


B ^ 4 A C 
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a source point exists. In the case of a quadratic in y, 


B =(D V Kd V 3 - 2(V /V )((V /V )r - r 
rk ry-' ' rs rvM' rz'^ rv^ v 


rs ry r rz ry y z 


A = 


C = 


1 +(V /V Y‘ 

\ rz ry/ 


D^(r /d)+ (V /V 

X ■ rz ry y 

2 2 

-D - (D r V /d V ) 

Y rx ry 
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It is now possible to outline the scheme for generating 
data with the computer model. 

a. ) the moon is advanced in the orbit by incrementing 
the time, thus determining the position r(t), and the deteC" 
tor look direction D Ct) . 

b. ) at each position, the values of t cal- 

culated based on B and V^. . 

c. ) the tests ,a\*h^cosS r and B>4 A C are per- 
formed and if they are both satisfied, dN ) is stored 

in the appropriate energy channel. 

To allow for the higher band pass at higher energies, 
all particles from energy^— to^4"^^ that satisfy the 
conditions described above, should be counted in the energy 
channel Due to the large amount of computer time requir- 

ed to test a nearly continuous range of energies, only the 
twenty center energies were used and each dN ( ^ , a ) was 
multiplied by 10% of the center energy. This is in keeping 
V7ith the qualitative nature of the model. 

Hidden in this technique is a numerical integration 
of the source points over the surface of the bow shock. That 
is, by making the above tests at each position r, over a wide 
range of energies and pitch angles , each source point that 
contributes a particle that satisfies the tests is counted. 
This can be seen by examining the special case shown in figure 


4-5b. 
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Chapter 5 

5.1 Inf roduc t 1 on 

Two basic approaches were followed in the analysis of 
particle data and interplanetary magnetic field data to deter*' 
mine whether or not the physical processes, discussed in 
chapter 3 are in operation in the pre and post bow shock re- 
gions of the lunar orbit. 

The first approach consisted of searching for the ex- 
pected correlation in interplanetary field configuration and 
observability of particles with the detector on the lunar 
surface while the second approach consisted of attempting 
to generate synthetic SIDE data, using the computer model 
and to d,iscovea; in vrhich ways these data resembled the real 
data. 

Simultaneous particle and magnetic field data were 
available from the Apollo 14 SIDE and Explorer 35 magnetome- 
ter respectively , for three periods of interest during May 
and June of 1972. Ths magnetometer data tapes v/ere provided 
by D.S. Colburn, Ames Research Center , Mof f et Field, Califor- 
nia . The regions in the lunar orbit corresponding to these 
periods and their associated times are shov/n in figure 5-1. 
Discussions of the analysis and results of the analysis of 
these data follow under the subtitles of Real Data and 


Synthetic Data. 
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Figure 5-1. The three regions of the lunar orbit 
in which magnetometer and ion data were collected 
for the analysis done in this thesis. The distance 
between the tick marks represents the distance 
traveled by the moon in 1 day. 
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5.2 Real Data 

In order to determine whether or not any correlation 
between B field configuration and particles observed exists 
in the data, it is necessary to devise a technique of com- 
paring the known B field components with the known particle 
accumulations of the SIDE for the same time interval. 

The obvious thing that comes to mind is to simply plot • 
the component of the field of interest and the integral par- 
ticle flux for the same time interval and visually inspect 
these plots for any correlation. An attempt at this proce- 
dure was made and the result is shown in figure 5-2. Al- 
though a weak visual correlation is evident from these data, 
it is obvious that no concrete conclusions can be drawn 
from this type of analysis. In order to overcome this problem 
we could divide this plot into small time intervals and con- 
sider each interval separately, noting the value of the field 
latitude and the number of particles observed at that lati- 
tude for each interval. That is, we would be interested in 
constructing from, the data, a plot of the average number of 
particles observed when the interplanetary B field has a 
particular configuration, vs, the field configuration. 

Such an analysis was done using the SDS 910 computer 
and data tapes from the Apollo 14 SIDE and the Explorer 35 
magnetometer. 
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Correlations in solar ecliptic latitude of the field 
and solar ecliptic longitude of the field, separately, were 
sought as were correlations with the field latitude and 
longitude considered jointly. Details of the computer analy- 
sis can be found in the appendix on computer programs. 

5.3 Results from the Real Data 

The results of the analysis using the real data are 
shown in figure 3 through 12 of this chapter, for the three 
periods of interest, which will henceforth be referred to as 
the inbound, tail and outbound regions. The inbound region 
corresponds to the pre-shock crossing region and the out- 
bound region corresponds to the post-shock crossing region 
of the lunar orbit. The inbound region consists of two days 
of data taken prior to entering the magnetosheath, beginning 
at 141 le^. The tail region begins at 148 0 and continues 

for two days? this period is entirely in the geomagnetic tail 
and the analysis in this region was done as a check on the 
validity of this technique. The outbound regions begins a- 
bout 3 hours after shock crossing at 153*^ 19^ 24^ and ex- 
tends for three days toward the sunset te3nninator. 

The format of the first three graphs is designed to show 
the average number of counts observed at a given solar eclip- 
tic latitude of the fie3.d vs, the latitude and the frequency 
of occurence of that latitude in units of, number of times 
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Figure 5 ”2 , Visual correlation of field latitude 
and integral flux, A correlation is especially 
apparent in the second half of this event. 

Figure 5-3 . Correlation of average integral counts 
and SE latitude of B over the inbound region. Bot- 
tom part of the graph shows the number of times a 
particular value of the latitude v;as encountered 
in the two day period over which this analysis was 
done . 

Figure 5-4. Lack of correlation in field latitude 
and counts in the tail region is apparent in this 
figure . 

Figure S-5, Correlation of field latitude and counts. 
Note the broad central peak indicating that particles 
are preferentially observed when the field is near 
the ecliptic plane. The enrichment of the high 
positive latitudes of this distribution is due to 
lunar atmospheric ions. 
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observad per number of days of observation vs, the latitude. 
That is the lower portion of each graph displays simply the 
ntimber of times a particular value of the field latitude was 
observed during the two or three day period of observation 
and the upper portion of each graph displays the average 
number of particles observed when at that latitude. Average 
integral counts will be taken to mean the average that is ob- 
tained by dividing the total number of particles accimulated 
at a particular latitude over the period of observation by 
the total number of times that that particular value of the 
latitude was observed during the same period. 

Due to the fact that the energy bandpass of the detector 
is proportional to the center energy, using counts instead 
of integral flux does not always provide a clear physical 
picture from the data. In view of this, the distributions 
of observed flux vs . field latitude v;ere calculated for the 
three regions as a check on the validity of using counts for 
this analysis. The results are shown in figure 5-6. By 
comparing these results with those in figures 5-3 through 
5-5, it is clear that using counts does not substantially 
change the form of the distributions. Comparison of the two 
results for the inbound region does show a fairly large in- 
flation of the high negative latitudes in the integral fltix 
representation that is not as apparent in the representation 
using counts. The reason for this is that the events that 
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caused the inflation have lower energies L ^ < 500 eV) than 
those of the broad central peak. Using integral flux tends 
to amplify the lower energies due to the smaller bandpass 
at these energies. The particles found in the upper and 
lower ends of the distributions of figures 5-3 to 5-6 are due 
largely to lianar atmospheric ions that are accelerated by 
the interplanetary electric field and the effect of these 
particles on the present analysis is discusseo more fully 
in chapter 6. 

The distributions shown in figures 5-7 to 5-9 are analo- 
gous to those in figures 5-3 to 5-5 where the solar ecliptic 
longitude of the field was the parameter of interest. It 
should be noted that the analysis was carried out here us- 
ing the uncorrected values of the field longitude. This 
has no bearing on the physical results , and in any physical 
interpretation of these results/ the scale in longitude may 
be simply shifted by 28°, i.e. the corrected value of the 
longitude is obtained by adding 28° to the uncorrected value. 
For a discussion of the error in longitude and the technique 
used to correct for it , see chapter 2 . 

The last three graphs of this series display the results 
of considering the latitude and longitude jointly. The eval- 
uation of the parameters in this treatment is over wider in- 
tervals, 10° latitude by 20° longitude, but the sajne proce- 
dure was used to determine the value in each two dimensional 
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Figure 5-6, Correlation in field latitude and int- 
egral flux. High latitude events are enhanced by the 
geometric factors, (see Appendix 3) . 

Figure 5-7. Correlation of field longitude with counts , 
for the inbound region. The distribution exhibits a 
double peak but is not well defined due to the scarcity 
of data for this region. 

Figure 5-8. Unexpected correlation in field longitude 
and counts in the geomagnetic tail. Note the well de- 
fined peaks in the frequency' distribution and the differ- 
ence in height of the two major peaks in the counts dist- 
ribution . i 

Figure 5-9 . Correlation in f ield longitude and counts 
for the outbound region. The two major peaks are about 
180° apart. 

Figure 5-XQ through 5-12. Correlations of field latitude 
and longitude considered simultaneously with average in- 
tegral counts. These plots contain basically the same 
information that is in figures 5-3 through 5-9. 


n 1 

INBOUND 

pT nr«IWi 



I 



-90 0 +90 


SE latitude of B 


Figure 5-6 







FREQUENCY 



Figure 5-8 


W INTEGRRL CTS 


:ncy 


\ 

\ 

! 



Ri/ INTEGRRL CTS 












94 


interval. The average integral counts in each interval of 
the grid is represented by the number of lines in that inter- 
val where each line represents five or ten counts. The ab- 
solute value of any of the numbers is of no consequence and 
no attempt is made here to display the numbers in a rigorous- 
ly quantitative fashion. These graphs were also produced 
with uncorrected longitudes and the same correction angle 
may be applied to them as was done for the case involving 
only the longitude correlation. 

After generally surveying these results, the first im- 
portant conclusion that can be drawn is that there is indeed 
a correlation between interplanetary field parameters and 
particle observations by the SIDE at the lunar surface . This 
statement can be made on the assumption that xf the number of 
particles observed by the detector did not depend on the 
configuration of the field, then the average number of parti- 
cles observed at any given value of a field parameter would 
not be very different from the average number of particles 
observed at any other value of the field parameter. In 
fact, that is the reason for carrying out the analysis in the 
tail region. The mechanism of the model is not expected to 
be operative in the tail region and a lack of correlation in 
field latitude can be seen in figure 5-4. Although the lati- 
tude of the geomagnetic tail field is sharply pealced about 
0° in frequency of occurence, the average number of particles 
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observed at any ef the latitudes is not too different from 
that observed at any other latitude. A clear lack of corre- 
lation in geomagnetic tail field longitude was not observed 
and this somewhat unexpected result will be discussed in a 
separate subsection below. 

With the foregoing in mind, a correlation of interplane- 
tary field latitude is apparent in figure 5-5 for the out- 
bound region. A rather broadly peaked distribution is ob- 
served with the peak centered around -10 to -15 degrees, 
i.e. with the interplanetary B field near the ecliptic 
plane. Although the peak is not sharp, for reasons that will 
be discussed later, there 'is a tendency for more particles 
to be seen, on the average, when the field latitude has a 
value that is in the neighborhood of 0”, i.e. with 


-25° < ®lat ~ +15°. Particles are still seen, however, at 
high positive values of field latitude and these are probably 
due to the outwardly directed interplanetary E field at high 
positive latitudes that accelerates the particles directly 
into the detector'. Csee fig, 3-4) 


The results for the inbound region, figure 5-3, a^e not 
as well defined, due primarily to the shorter period of ob- 
servation, but again there is a broad peak that centers near 
0° latitude with some particles still being seen at high neg' 
ative values of the latitude. In the inbound region, if the 


field latitude has high negative values, the interplanetary 
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"e field again points out into the detector on this side of the 
raagnetosphere Csee figure 3-4) accounting for these high neg- 
ative latitude particles . 

At this point it is possible to conclude that the ob- 
servability of particles in the inbound and outbound regions 
of the lunar orbit is dependent on interplanetary field lat- 
itude and the dependence does roughly correspond to that ex- 
pected from the model. In addition, there is indication from 
the data that the asymmetric interplanetary E field at high 
positive and negative values of the B field latitude is pro- 
viding some of the particles observed by the detector in 
these regions. 

Figure 5-9 shows a correlation in field longitude for 
the outbound region with a fairly obvious doubly peaked dis- 
tribution, the first peak coming at about 84“ uncorrected 
and the second smaller peak coming at about 264° uncorrected, 
a separation of 180° in longitude. One would also expect 
from the m.odel a doubly peaked distribution with approximate- 
ly the same size peaks, separated by about 180° since the 
polarity of the field is unimportant in determining the re- 
sultant guiding center trajectory of a particle for any par- 
ticular value of the field latitude. For example, if the 
field latitude is say -10° while the longitude is 84° and 
+10° V7hen the longitude is 264°, the same physical situation 
exists with respect to the model and the polarity is unim- 
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portant. If, hoxvever, the latitude is also -10® at the 264“ 
position in longitude, a different physical situation doe 
exist and the guiding center trajectories will depend on the 
polarity of the field. Figure 5-13 illustrates this effect. 
Inspection of figure 5—12 reveals that the latitudes of the 
peaks in figure 5-9 are not too different. That is, differ- 
ent physical situations exist for the two peaks. The plane 
of the guiding center trajectories of case (1) in figure 
5-13 makes a smaller angle with the plane of the lunar orbit 
than does the plane of guiding center trajectories of case (3), 
possibly accounting for t> i favorability of observing parti- 
cles in cased) , i.e. the higher peak at about 84® longitude. 

The data from the inbound region, figure 5-7, also show 
correlation in longitude of the field with the number of par- 
ticles observed, but again, the shorter period of observation 
results in a less well defined distribution. Despite this, 
the distribution also shows a double peak with a less pro- 
nounced difference in height. A less pronounced difference in 
latitude for the two peaks is found by examining figure 5-10. 
The two peaks occur at latitudes much nearer the ecliptic 
plane thus lessening the effect shown in figure 5-13. These 
differences are consistent with the ideas discussed above 
concerning the actual physical situation occuring for the 


two different peaks . 
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Figure 5-13. Views in two different planes of the phy- 
sically different situations ocurring for field config- 
urations in which the longitudes are 180° out of phase 
and the latitudes are equal. 



Cases 1 and 2 physically equivalent with respect to model 
Cases 1 and 3 not physically equivalent 

Cases 1 and 3 are favorable to observation of particles 
but case 1 is more favorable than case 3 according to the 
data in figure 5-8. 


Figure 5-13 
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The results of the longitude analysis show that there 
is correlation in observability of particles with the SIDE 
and the solar ecliptic longitude of the interplanetary B 
field and this correlation further supports the model dis- 
cussed in chapter 3 . 

5 . 4 Interpretation of tail passage 

As mentioned above, figure 5-8, eKhibits a doubly peak- 
ed distribution with the peak corresponding to the field of 
the geomagnetic tail pointing in an anti-solar direction at 
about a factor of 2 greater than the peak corresponding to 
the tail field pointing in the solar direction. This may not 
be significant, but since there is no obvious lack of correla- 
tion in longitude and number of particles observed, the phy- 
sical situation in this region should be examined more closely. 

At the beginning of the tail region, (148^ 0^ l'^) , the 
solar ecliptic longitude of the geomagnetic tail field is 
fairly stable around 358°, corrected, or with the tail field 
approximately in the solar direction. At approximately 
148°^ 5^, the field longitude changes abruptly and becomes 
stable about the value 178°, or the anti-solar direction and 
remains there throughout the two day period. The fact that 
no more abrupt changes in field longitude, i.e. nautical sheet 
crossings, are observed over this period can be understood by 
looking at a diagram of how the relative distances from the 
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ecliptic plane of the moon and the neutral sheet at the long- 
itude of the moon vary with time over the tail region at 

60 R , The curve R in figure 5-14 represents the distance 
0 

of the neutral plane from the ecliptic plane at the longitude 
of the moon, taking into account the seasonal offset and the 
diurnal variation; an off axis pivot point is also assiimed 
due to the effect of the solar wind aberration angle of •^4'’. 
The curve d of figure 5-14 is simply the distance of the moon 
from the ecliptic plane at 60 R as it moves toward higher 
longitudes. The sense of the geomagnetic tail field, in the 
region in figure 5-14 that is above the curve P., is in the 
solar direction and the field has the opposite sense in the 
region beneath the curve R. The motions in this diagram 
should be considered from the point of view of an observer 
on the earth, looking in the anti-solar direction. The moon 
spends a short time in the northern lobe region of the tail 
field and then crosses into the southern lobe for the rest 
of this period. 

All that can actually be concluded from the result in 
figure 5-8, is that there is a tendency for more particles to 
be seen by the detector on the average, when the tail field 
points in the anti-solar direction, i.e. toward the detector, 
than when the field points in the solar direction, away from 
the detector. Note that the detector looks earthward in this 
region of the orbit. There is evidence here for a shadowing 
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Figure 5-14 . View from earth, looking in the anti-solar 
direction. Cuirve R represents the distance in of the 
magnetic neutral plane from the ecliptic plane as a func- 
tion of time. Curve d represents the path of the moon 
with respect to the ecliptic plane for the same time per- 
iod. Note that the moon spends most of the two day period 


in the southern lobe of the earth's field. 
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effect by the moon. Van Allen et al. I197U , have suggest- 
ed a model by which thermal anisotropies in the solar wind 
can provide an asymmetric access of protons to the magneto- 
sphere, The basic idea of this model is that interplanetary 
field lines connect to geomagnetic field lines in the dis- 
tant magnetotail (~100 . When the interplanetary field 

has a southward component and is directed av/ay from the sun 
it connects to the tail field lines as shown in figure 5-15. 

It can be seen from examination of this figure that the pos- 
sible shadowing seen in the present analysis is consistent 
with the Van Allen model. By examining other features of 
the data in this region, for example the enei-gy spectra, a 
better determination of the origin of this phenomenon could 
be made. 

5 . 5 Synthetic Data 

The basic philosophy underlying the second approach has 
been described in chapters 3 and 4 and the results of the ap- 
plication of that philosophy and a brief restatement of the 
scheme for generating the synthetic data are given in the fol- 
lowing sections . 

5.6 Brief restatement of the scheme 

a.) the position vector ^(t) and the detector look dir- 
ection D (t) are calculated. 
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Figure 5-15. Van Allen model for asymmetric access of 
solar wind protons into the magnetosphere. The clustered 
arrows on the field lines represent proton thermal an- 
isotropies in the solar wind. [Van Allen et al. (1971)]. 





104 


b, l at each value of t, the resultant velocity vectors 
^ f o. ) are calculated based on B and 

c. ) two constraints are then imposed to consider the par- 
ticle countable. 

1. ) the angle between V^( ^ / u) and D is required to 
be less than some limiting angle, the value of which is the 

subject of the next subsection. 

2 . ) the particle must have a source point on the 

shock surface. 

5.7 Detector Looi< Cone 

The detector look cone will be taken to mean the solid 
angle centered on the detector look direction D and having 
half angle 6 . The nominal value of 6 is as deter- 

mined from the angular calibration of the instrument using 
a parallel ion beam in the lab. It must be kept in mind that 
when the detector is observing ions of large gyroradius com- 
pared to the dimensions of the entrance aperture of the in- 
strument, the look cone of 3° half angle does not prevent the 
detector from seeing particles whose guiding center traject- 
ories come from directions that are not in the field of view 
of the look cone. This is especially important in interpret- 
ing the data in this chapter since the gyroradii of bow 
shock protons are typically hundreds of kilometers and the 
guiding center trajectory does not represent the actual path 


105 


of the particle. One would expect that if the SIDE could see 
particles only in the 6° look cone that the distributions 
shown in figures 5-3 through 5-8 would have sharp peaks about 
6® wide. One would also expect that if the effective size 
of the look cone was say 8 ~ 90°, then when the detector was 
very near the shock it would always see particles since the 
angle subtended at the detector by the shock surface would be 
very large. The distributions are very broadly peaked in- 
dicating that the effective cone angle for observing bow 
shock protons is much larger than the nominal 6 , but that 
the half angle is not as large as 90°. 

The influence of this effectively larger cone angle on 
the computer model is to relax the criterion that the resul- 
tant velocity vector lie within a 3° half angle of the detec- 
tor look direction. It was determined from testing the com- 
puter model that a cone half angle of about 60° gave the best 
agreeipemt with the real data. For further discussion of the 
detector look cone and its relation to a particle ' s pitch 
angle see the discussion of the range of observable pitch 
angles later in this chapter. 

5 . 8 Application of Computer model with idealized conditions. 

The first application of the computer model was done 
using the Rice IBM 370 computer, using artificial input para- 
meters, with the exception of those parameters relating to the 
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lunar orbit. Explorer 35 data was not available at this 
time so typical values of the interplanetary field parameters 
and solar wind velocities were used to test the qualitative 
behavior of the model. The result, shown in figure 5-16, was 
generated using a solar wind velocity of 400 km/sec, a field 
latitude of 0® and a field longitude of 315®, Ninety pitch 
angles were tested at each energy and a spectrum was generated 
for every 40 minutes of lunar travel for the 5 day period 
shown . 

Particles with energies of less than about 500 eV were 
not observed and only the highest three or so energy channels 
(3 to 3,5 keV/q) contained counts when the last spectrum was 
generated. The cutoff shoulder moved steadily toward higher 
energies with increasing time for each spectrum some particles 
were seen. The behavior of the model exhibited in figure 5-16 
is roughly characteristic of what is actually observed dur- 
ing the post bow shock crossing regions one each lunation. 

The spectra from the real data just after bow shock crossing 
are generally peaked around 500 eV to 1250 eV and after 4 or 
5 days only higher energy particles are seen if any at all. 
Figure 5-17 shows a typical post bow shock region for compari- 
son. 

It should be noted that in this highly idealized case, 
at least the gross features found in the real data are re- 
produced. The low energy peaks in the real data are due to 
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Ficnire 5-16 . Synthetic data using typical interplanetary 
field conditions, namely ®long~^^^°* 

energy channels of the instrument are plotted on the hor- 
izontal scale representing 10 eV to 3500 eV; log^Q of the 
average counts is plotted vertically and time runs up the 
oblique axis. 
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Figure 5-16 


Figure 5-17 . 


Typical SIDE counting rate spectra from the 
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lunar atmospheric ions. 

The solar wind velocity and the interplanetary magnetic 
field are highly time varying quantities and only if this is 
taken into account in generating these spectra can one expect 
to effectively test the model of chapter 3 by generating 
synthetic data, A more elaborate procedure for generating 
synthetic data would be to calculate actual particle traject- 
ories which may be done in future analyses of this problem. 


5*9 Inclusion of interplanetary magnetic field data. 

When Explorer 35 magnetometer data became available for 
the outbound region^ dxscussed above ^ the computer model 
was refitted for the SDS 910. Epheineris data tapes for 
Apollo 14, containing information on lunar position and de- 
tector look direction also became available so that it was 
no longer necessary to calculate the orbital parameters as 
described in chapter 4. The field information from the mag- 
netometer and the position information from the ephemeris 
were synchronized and fed into the model and analyzed at 8 
minute intervals using the scheme discussed above. The solar 
wind velocity was still taken to be 400 km/sec throughout. 

The counting rate spectra from the computer model, at 
8 minute intervals, for the outbound region are shown in 
figures 5-18 and 5-19, The gaps in time are due to the non- 
continuous data from the Explorer 35 magnetometer. It is 


Figure 5-18. Synthetic data generated using the real 
values of the magnetic field parameters and the real 
ephemeris data as input to the program (see Appendix 3) . 
Data gaps are due to the gaps in the Explorer 35 tapes. 

Figure 5-19. Continuation in time of figure 5-18. 

Figure 5-20. Real ion counting rate spectra for the 
time corresponding to that in figure 5-18, The high 
sharp peaks are transmission noise. 

Figure 5-21. Continuation in time of figure 5-20. 




Figure 5-18 


APOLLO 14- 12 

START TIHE 155 
6 ni:S IMTl/^LS 
12 . 0«3HR/Oll^ 



1 

EMFRGY (eV) 


Figure 5-19 


APOLLO l«l- 12 

START TIME 153 19 25 

8 MIM 

IE . 00HR/DII^ 


5. 


APOLLO 14- 12 

START TIME 155 3 
8 MIN Al<^ 

12 .00HR/DIU' 



1008 

ENERGY (eV) 


3 5 > 1 1 ' 


Figure 5-21 


Ill 


inmiediately obvious that the inclusion of the magnetic field 
data causes more structured synthetic counting rate spectra, 
but again it can be seen from examining these results that 
the observed particle energies are in the higher regions of 
the spectrum. That is, particles with energies less than 
about 500 eV are generally not observed. It is also interest- 
ing to note that the spectra are no longer always continuous 
in energy after reaching the cut-on energy; multiply peaked 
spectra are readily apparent. The results shown in figures 
5-18 and 5-19 resemble real counting rate spectra much more 
closely than those of figure 5-16, even without including 
the changes in solar wind velocity. 

The real counting rate spectra for the outbound region 
are shown in figure 5-20 and 5-21 for comparison. Although 
there are some obvious differences, the data generated by the 
computer model are not completely unlike the real data with 
both exhibiting structured spectra with counts mostly in 
the higher energy channels. The noisy spectra observed in 
the early post bow shock region of this lunation is not al- 
ways present and could be due to some turbulent effect, seen 
only while the moon is close to the shock itself or possible 
temperature effects on the instrument are responsible for 
this noise. Unfortunately, event by event comparisons can- 
not readily be made due to the time gaps. 
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After generating the data shown in figures 5-18 and 
5-19, an analysis, using the real data, was applied and this 
result is shown in figure 5—22 . Due to the fact that the 
spectra were generated only every 8 minutes, the n\unber of 
times a particular value of the field latitude was observed 
is considerably less than the number observed in the analysis 
of the real data and hence, the distribution of average in- 
tegral counts is not as well defined. The results for the real 
data for the same period are shown in figure 5-23 for compar- 
ison. 

It may seem circular to analyze the generated data under 
the philosophy for the first approach since the model was 
forced to operate in generating these data whereas in analyz- 
ing the real data, we were looking for evidence that the mo- 
del was operating. It is, nevertheless, interesting to note 
that by forcing the model, a distribution is obtained that 
is very much like the actual distribution. 

5.10 Range of observable pitch angles. 

Due to the rather narrow field of view of the detector, 

(3® half angle) , and the finite extent of the bow shock sur- 
face, some particles leaving the shock surface with large 
pitch angles may not be seen by the detector. It is possible 
to estimate the range of observable pitch angles at any given 
position in orbit by considering the simple geometry of fig- 


Figure 5-22. Correlation in field latitude and counts 
using synthetic ion data, as input. 

Figure 5-23. Correlation of field latitude and counts 
for the same time period using real ion data. 
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ure 5-24. The actual bow shock surface is approximated by a 
flat surface subtending approximately the same angle at the 
detector that the actual surface would subtend. The points 
L represent the limits of the source region, D represents 
the radius of a cross section of the bow shock in a dawn- 
dusk meridianal cut and R is the geocentric distance of the 
moon. S is the half angle of the detector look cone, y is 
the angle subtended at the moon by the bow shock and a 's 
represent pitch angles. 

It can be seen from the figure that a particle leaving 
the shock at point C must have a pitch angle of S or less to 
be detected by the instrument on the lunar surface. However, 
if the particle comes from the point L it must have a pitch 
angle in the range 

y-8<a<y + 8 

to be detected. Although each point on the surface defines 
such a range of pitch angles, the particles coming from the 
points L define the upper limit of the pitch angle range 
for a given distance from the shock. That is, particles 
with pitch angles 

a > y + 8 

cannot get into the detector from any point on the shock at 
the distance defining y . 


Figure 5-24. Illustration of the relationship between 
the detector look cone angle, the angle subtended at the 
moon by the bow shock and the pitch angle of the particle. 
The figure shows the moon at last quarter in the dusk-dawn 
meridional plane. The dotted particle trajectory represents 
the upper limit to the pitch angle at this distance. 
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When the moon is very near the shock y is very large 
and the particles of high pitch angle may be observed while 
as the moon moves away from the shock y gets smaller, limit- 
ing the range of observable pitch angles. The limit of ob- 
servable pitch angles for a case in which the a -range is 
fairly restricted can be estimated by calculating y when the 
moon is on the dusk-dawn meridian. For this case, 

D ^ 20 R 

e 

R 60 R 

e 



a <23° 


This limit is probably one of the more restrictive due 
to the fact that as the moon moves further toward lunar 
nightime the detector sees a larger shock surface, thereby 
increasing y and as mentioned above, y is larger when the 
moon is nearer the shock surface . 

The effect of this limited range of etch angles was not 
considered in constructing the computer model, but the effect 
of this neglect should not be significant. Particles 
with high pitch angles have small // velocities causing their 
resultant guiding center trajectories to be more in an anti- 
solar direction than ones with small pitch angles indicating 
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that as one gets further away from the shock it becomes more 
difficult to see a high pitch angle particle based on consid- 
erations of guiding center trajectories alone. Rot reject- 
ing a particle in the model because of its high pitch angle 
does not significantly degrade the synthetic data since the 
particle would probably not be seen anyway because of its 
largely anti— solar guiding center trajectory. 

The above considerations neglect the fact that the 
pitch angle is not constant with regard to the guiding center 
trajectory, but similar behavior of the pitch angle limits 
can be expected if this effect is taken into consideration. 

5.11 Energy density of upstreaming protons 

To estimate the effect of the upstreaming protons on the 
incident solar wind it is appropriate to calculate the energy 
density of the particles and compare this with the energy 
density of the incident solar wind protons. 

By measuring the counterstreaming flux of the particles 
seen at the detector and assuming that each observed particle’s 
energy is due to its linear motion into the detector, one 
obtains an energy density of 

-13 3 

~ 7 X 10 ergs/cm 

S • s * 

for the 1 to 3.5 keV particles seen on a typical davm side 
bow shock region. The incident solar wind energy density is 
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—9 3 

of the order of 7 X 10 ergs/cm which implies that 

1? / 7 } ~ 10“^ (1 kev to 3.5 kev) 

B.S. S.W. 

One may conclude that the bow shock protons observed by 
SIDE would probably not cause a significant perturbation of 
the incident solar wind. 

Other investigators have calculated energy density for 
the counterstreaming protons that are considerably higher 
than this value, however none of these other values are 
greater than about 1% of the incident solar wind energy 
density for quiet conditions. Lin and co-workers [1974] 
calculate a value of 10 ergs/cm for protons in the range 
29 to 100 keV. Evidence from their data indicates that there 
is a sharp cutoff in counters treaming protons above 100 keV 
which implies that 100 keV is probably the upper limit in 
energy for bow shock protons. Asbridge et al. [1968] calcu- 
late an energy density of 10 ergs/cm in the range of 3 

to 6 keV which is 1% of the incident solar wind energy 
density. Scarf et al. [1970] have observed values of the 
counterstreaming energy density that are comparable to the 
incident solar wind energy density at times. However, it 
seems that large scale, long time perturbations of the inci- 
dent solar wind by the counterstreaming protons is not a 
significant factor in determining the general behavior of 
these particles. 
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5.12 Conclusions . 

The results of the analyses presented in this chapter 
give rise to the follovjing conclusions; 

a. ) the observability of particles at the lunar surface, 
during the pre and post shock regions of the lunar orbit, is 
at least partially dependent on the configuration of the in- 
terplanetary magnetic field. 

b. ) the evidence from the analysis of the data and 
from the behavior of the model suggests that the model dis- 
cussed in chapter 3 is probably operative in the pre and 
post shock regions of the lunar orbit and is responsible for 
the moderate to high energy particles observed in these re- 
gions. This result is not as well established in the pre 
shock region as in the post shock region and more data are 
needed to further analyze this region. 

c . ) The non-zero gyroradius of the particles cannot be 
ignored in the computer model but the effect it has on gen- 
erating synthetic data can be approximated by enlarging the 
detector look cone angle to include a larger range of guxd- 
ing center trajectories, 

d. ) the energy density of the counter streaming protons 
is not sufficient to significantly perturb incident solar 


wind conditions . 
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Chapter 6 


6 . 1 Disctission 

There are possible sources other than the bow shock of 
moderate to higher energy protons, ^ > 750 eV, seen in the bow 
shock regions, viz. 

a. ) lunar atmostpheric ions 

1. ) interplanetary E field ions 

2. ) surface E field ions 

b. ) solar wind ions and suprathermal solar wind ions 

c. ) cosmic rays 

d. ) unknown sources such as lunar bow shock protons 
The "bow shock region" will refer to the regions of the 
lunar orbit extending from shock crossing to about a day 
or so past the sunset terminator on the dawn side and the 
analogous region before shock crossing on the dusk side of 
the magnetosphere . 

The reasons that it is unlikely that these sources can 
account for the observations completely will be pointed out 
below. 

a.) lunar atmospheric ions-As discussed by Lindeinan [19 73] 
and Manka [1972] , ions from the lunar atmosphere can be ac- 
celerated into the detector by the interplanetary electric 
field when the field lies in the ecliptic plane. By 
considerations of surface density and temperature, ion scale 
height, ionization rate, and strength of the electric field. 
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the expected energy spectrum observed by the SIDE is given by, 

F( ^ [" 6-1 

where, 

n =surface number density 
p=ionization rate 
e=charge on the ion 
H=scale height 



^ =ion energy 

Using appropriate values of these parameters for the 
Apollo 15 event observed on day 312, 1971 Lindeman [1973] 
has determined the predicted flux values and compared them 
to the observed fluxes. The result is shown in figure 6-1 
(after Lindeman [1973]). The flux drops off rapidly after 
about 250 eV indicating that the higher energy ions observed 
in this same region are not lunar atmospheric ions which 
have been accelerated by the interplanetary electric field. 
It should also be pointed out that the distribution obtain- 
ed in figure 5-5 indicates that most of the particles were 
observed when the interplanetary electric field was not in 
the ecliptic plane in contrast to the results obtained by 
Lindeman [1973] , (figure 6-2) , and Manka [1972] , (figure 6-3) , 
that lunar atmospheric ions are seen when 


is in the 
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Figure 6-1. Predicted and observed energy spectra of 
lunar atmospheric ions. Note the sharp drop in energy 
after about 250 eV. [after Lindeman (1973)]. 

Figure 6-2. Integral flux vs. angle out of the ecliptic 
plane of the interplanetary *E field. The integral flux 
drops smoothly as the field goes out of the ecliptic. The 
bottom part of the graph displays ^ /1 f| vs. angle out 
of the ecliptic plane. [Lindeman, 1973] 

Figure 6-3, Correlation in integral flux and ® 

implies an E field that is about in the ecliptic plane 
and pointing into the detector. [Manka (1972)] 
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ecliptic plane. 

Indications are that at the terminator the largest value of 
an attractive lunar surface potential is on the order of -100 V 
[Fenner et al,. , 1973 ] , which implies that essentially thermal 
atmospheric ions could not attain energies much greater tlian 
about 100 eV before being detected by the SIDE. The lunar 
surface potential is obviously not responsible for the part- 
icles of energy greater than about lOO eV. It is difficult 
to imagine any other local electric fields of sufficient 
strength to accelerate atmospheric ions to the observed ener- 
gies . 

b.) solar wind ions-The thermal energy of a solar wind 
proton is on the order of 10 ev and its energy due to bulk 
flow motion, {V ^ 400 km/sec) , is > 1000 eV. The angle 
between the look direction of the detector and the solar v?ind 
flow in the bow shock region eventually becomes >9 0°, yet the 

high energy flows are still observed. 

Although local lunar surface magnetic fields as 

high as 300 gammas have been measured, Dyal et al. [1972] , 
these fields are not of sufficient strength to deviated solar 
wind particles into the detector. In addition, there is still 
no known acceleration mechanism that could produce the 3500 eV 
particles that are observed. There is no reason to suspect 
that the spectrum ends abruptly at 3500 eV as well. Counter- 
streaming protons of energy 3 to 6 keV have been observed in 
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this region of space by other investigators, viz. Asbridge 

et al. [1968]. Lin et al., [1974] have observed counter- 

streaming protons in the energy range 30 to 100 keV . Although 
solar wind protons in the energy range 5 < ^ < 50 keV known 

as suprathermal solar wind protons have been observed [Frank, 

1970] , these energetic proton events are not characteristic of 

the quiet solar wind. Since these suprathermal protons come 

from the solar direction one would again have to invoke some 

mechanism to deviate the particles into the detector . 

c.} cosmic rays-Very high energy protons ( ^ > 50 MeV) 

can penetrate the case of the SIDE instrument and activate 
the channeltron causing counts to appear in whatever energy 
channel is being sampled at the time. The cosmic ray flux 
in the interplanetary medium is responsible for part of the 
background noise during instrument daytime and is responsible 
for most of the background during instrument nightime . These 
counts in random channels are generally very low (1 or 2) , and 
are only enhanced during unusual solar proton events such as 
solar flares as discussed by Medrano [1973]. That is, cosmic 
ray protons produce a random background very low intensity 
in quiet times and ir disturbed times the background counts 
rise independent of energy channel making their identification 


easy . 
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d.) other sources-At the time of this writing no lunar 
bow shock capable of producing 3 keV protons has been observ- 
ed. Other complicated solar wind-lunar interactions may occur 
to produce these particles, but such processes have not been 
reported . 

The above considerations alone indicate strongly that the 
high energy particles observed in this region are indeed 
protons from the bow shock since it is the most likely other 
source. Other corroborating evidence that the high energy 
particles are from the bow shock can be summarized as follows. 

a. ) Protons of the same energies as the observed ones have 

been previously observed in this region by several investiga- 
tors so it is not unreasonable to expect that the SIDE would 
detect them. The instrument looks back at the shock in this 
region of the orbit. 

b. ) Results of the data analysis of this thesis indicate 
that these high energy particles are influenced by the inter- 
planetary magnetic field in such a way as to be consistent 
with the proposed EXB drift model. 

It is apparent that bow shock particles are a dominant 
feature of the data observed in the bow shock regions which 
suggests a reclassification of the recurring ion events 
discussed by Lindeman [1971,1973]. Lindeman classified the 
ion spectra in the extra-magnetospheric regions of the lunar 
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orbit into five basic types, two of which are clearly ident- 
ifiable with bow shock protons. The Type I spectrum has a 
single peak from 2000 eV to 3500 eV and Type II has a double 
peak, the first peak from 30 eV to 250 eV and the second 
from 2000 eV to 3500 eV, (see figure 6-4) . Subsequent inves- 
tigation by Lindeman [1973], has shown that the low energy e- 
vents in the bow shock region can be classified into three 
types . 

a. ) VXB events 10 eV to 500 eV 

b. ) Lunar surface potenial (LSP) events 10 eV to 100 eV 

c. ) Instrument resonance events (IRE's) ^50 eV 

In light of these findings, the low energy peak of the 

Type II event is easily identified as being due to LSP and 
VXB events near the terminator with the high energy peak be- 
ing explained by the presence of bow shock protons in this 
region. The Type II event is merely a Type I event with the 
contribution to the spectrum of the LSP and VXB events . 

The ion data observed in the bow shock region can there- 
fore be classified as consisting of four basic types of events; 
a.,b.and c. above and 

d. ) Bow shock proton (BSP) events ^ a 750 eV 

Based on the present research it is possible then to state 
a general criterion for identifying BSP events in the data. 
During passage of the moon through the bow shock region, any 
observed particles having energy greater than or about equal 


Figure 6-4. Recurring ion events Types I and II . Type I 
has only a high energy peak v/hile Type II has a low energy 
peak in addition to the high energy peak, [redrawn after 
Lindeman (1971) ] 
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to 750 eV are probably bow shock protons / with the exception 
of clearly identifiable high energy events associated with 
unusual activity in the interplanetary medium. 

Figure 6-5 shows representative counting rate spectra in 
the bow shock region near the sunset terminator with two of 
the above types of events easily identifiable. LSP's and BSP's 
are readily apparent and a tentative identification of V X B 
events can be made . The marginality in the identification 
of the V X "B events could be easily removed by looking at the 
interplanetary field behavior for this period and remembering 
that B out of the ecliptic implies V X B events and B in the 
ecliptic implies BSP events. 
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Figure 6-5. Representative part of the bow shock region 
showing varied ion spectra near the sunset terminator. 
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6.2 Further Research 

Limitations on the accuracy of the present analysis are 
due primarily to a lack of sufficient simultaneous magnetic 
field and ion data. An analysis of the bow shock regions 
carried out over several lunations would much enhance the 
significance of the results and would likely improve the agree 
ment of the data with what is expected from the model. This 
can be seen by noting the better correlation achieved in the 
outbound region, (see figure 5-5) , where three days of data 
were analyzed as opposed to two days on the inbound side . 

The amount of simultaneous particle and field data available 
in any one lunation is also limited by the data gaps in the 
field data incurred while Explorer 35 is behind the moon. 

The lack of an onboard recorder causes data gaps of about 5 
hours on each orbit cf the satellite. 

Other problems exist such as the unknown spatial distri- 
bution of the field, an important factor in determining the 
resultant guiding center trajectory of a particle. Not much 
can be done about this problem but it appears that the var- 
iation of the field over a distance comparable to the proton 
gyroradius is not significant since particles do seem to make 
it to the moon via the EXB mechanism. Likewise the problem 
of including the effects of the local surface fields is very 
difficult in that the spatial configuration and extent is not 
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well known, but again these fields dp not seem to play a 
significant role from the self consistent point of view. That 
is, if one assumes a model and observes the behavior predict- 
ed by the model then other influences must not be too signifi- 
cant . 

In some instances the field is also highly time varying 
and no account of this has been taken in the present analysis. 
If the field is varying rapidly with respect to a character- 
istic time of the motion, say the gyroperiod, ( 12 sec) , the 

guiding center treatment breaks down. Thera were several in- 
stances in the data where the field was stable with respect 
to a gyroperiod but the inclusion of times where the field 
was varying rapidly would detract from the accuracy of the 
results. An obvious way to improve the results would be to 
confine the analysis to times of stable field conditions . 

As noted above, the bow shock reyr'ons are of a highly 
diverse character in ion spectra v/ith bow shock protons con- 
tributing a large portion of the observed ions but certainly 
not all of them. This would imply that any correlation of 
particles observed with field configuration would be further 
obscured by including particles whose observability did not 
depend on the field configuration. All events were included 
in this analysis. Obviously, the contribution by the LSP's, 
the ire's and any other event associated with the lunar atmos- 
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phere would be detrimental to the results. This point can 
be clearly seen by noting the effect of including all of 
the ion events on the distributions shown in figures 5-5 and 
5-3. la particular, the inclusion of the VXB events inflates 
the upper end of the distribution in the outbound region and 
the lov/er end of the distribution in the inbound region. For 
this particular type of event the contributions to the dist- 
ributions cannot really be considered detrimental since VXB 
events do depend on the field configuration and are easily 
distinguishable from BSP events. 

Again by considering isolated events where the field is 
stable and where particles from other known sources can be 
removed from the analysis, e.g. near the terminator, the cor- 
relation would be greatly improved. 

With the availability of good data from all three SIDE'S 
and Explorer 35, the directional nature of the proposed model 
could be more effectively tested. The wide variation in long- 
itude of the look direction of the three detectors would al- 
low a cross correlation in dependence on longitude of the ob- 
servability of ions. The use of three instriiments in this 
manner could also be helpful in determining the effective 
size of the look cone as discussed in chapter 5. 

Any investigation such as the one done for this thesis 

yields many new and interesting areas for further 
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work. Tv/o of the related areas of special interest are, the 
possible acceleration mechanisms associated with the bow shock 
itself and a study of the nature of the source points on the 
shock surface. If a very good fit to the real data can be 
achieved using a particular type of source function, then one 
is led to suspect that that source function probably approxi- 
mates the real one. Knowing the source function at the shock 
surface could be useful in determining which if any of various 
models for the acceleration mechanism is correct. Two of the 
major models for the acceleration mechanism are reviewed in 
the Appendices . Specular reflection at the shock and subse- 
quent energization by the interplanetary electric field has 
been proposed by Sonnerup [1969] and the mechanism of Fermi 
acceleration has been developed by Jokipii [1966] . 

However, the primary mission of the SIDE is to gather 
data that can be used to learn about the lunar atmosphere 
which suggests that the real value of this research and any 
further research in this area is that it aids in understand- 
ing the complicated spectrum of events occurring in the bov^ 
shock regions of the lunar orbit. These bow shock regions 
are the most fertile in the data for providing information 
on the lunar atmosphere and certainly warrant further inves- 
tigation. 
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Chapter 7 

7.1 Sungnary 

The basie objectives of the research done for this thesis 
were outlined in chapter 1 and will be briefly restated here. 

a. ) to determine if the particles of energy >750 eV ob- 
served in the pre and post bow shock regions are indeed bow 
shock particles. 

b. ) to establish a criterion for the identification of 
these particles. 

c. ) to propose and develop a physical model to explain 
the presence of the bow shock protons if they are found in 
the data . 

Such a model based on EXE drift trajectories was developed 
in chapter 3 leading to the following basic conclusions about 
the behavior of bow shock protons in this region. 

a. ) the observability of bow shock particles at the moon 
by this mechanism is highly dependent on the configuration 
of the interplanetary magnetic field. 

b. ) there is an inherent dusk-dawu asymmetry in the 
model that should be observable in the data. 

It was also pointed out that the present model/ although 
similar to that proposed by Asbridge et al. [1968] , is a more 
realistic one in that it describes all particles leaving the 
shock;,/ not just those with a special velocity and pitch angle 
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The details of the computer model and its connection to 
the physical model were discussed in chapter 4 v?here criteria 
for the observability of bow shock particles at the moon were 
established. For a particle to be observed at the moon by the 

-A -A 

EXB mechanism, it must, 

a. ) come from the right direction and 

b. ) have a source point on the shock surface. 

The plausibility of the physical model of chapter 3 was 
established by the data analysis, the results of which are 
displayed in chapter 5. The major conclusions drawn from 
interpretation of the data are, 

a. ) analysis, independent of the model, shows that the 
data behave as one would expect from the model . 

b . ) synthetic data generated by the computer model agree 
semi -quantitatively with the observed data. 

c. ) behavior of the computer model coupled with consider- 
ations of actual particle trajectories indicates that the 
effective size of the look cone of the detector is substanti- 
ally larger than the nominal value when applied to bow shock 
protons. 

Arguments against other possible sources for the observ- 
ed ions of energy > 750 eV are presented in chapter 6. 

These other sources are, 

a. ) lunar atmospheric ions 
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b. ) solar wind ions 

c . ) cosmic rays 

d. ) lunar bov; shock or other unknown sources 

and were rejected primarily on energy considerations. Lunar 
atmospheric ions and solar wind ions have much lower ensr~ 
gies in general than the observed moderate to high energy 
ions and it was pointed out that no known acceleration mech- 
anism at the moon is capable of producing particles of the 
observed energies. Cosmic rays were dismissed based on the 
well established character of penetrating particle events 
which are unlike bow shock proton spectra. 

The value of isolating and identifying bow shock protons 
in the data was also discussed in chapter 6 where it was noted 
that the usefulness of the bow shock region for studying the 
lunar atmosphere can be enhanced by separating BSP events fi:om 
the lunar atmosphere events. A general criterion for ident- 
ifying bow shock protons was given; particles of energy great- 
er than or about equal to 750 eV and found in the pre and 
post bow shock crossing are probably associated with the 
bow shock. 

7.2 . Conclusions 

a. ) The protons observed in the bow shock regions with 
energies > 750 ev are from the bow shock and can readily 
be identified in the data. 

b. ) The proposed EXB drift model is a viable one as long 
as its limitations are kept in mind. 
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Appendix 1 

Al.l Define an orthogonal coordinate system S and an orth- 
ogonal coordinate system S‘ that has the same origin as S 
but is oriented arbitrarily with respect to S otherwise. It 
is possible to make the axes of system S coincident with 
system S' by three sucessive counterclockwise rotations of 
the S axes. The angles through which one must rotate to per- 
form such a "transformation" of the coordinate axes are known 
as Euler angles. The angles , ijj , and 0 of figure Al.l are 
the Euler angles for the transformation from system S to S 
and the transformation matrix and its inverse are shown on 
the next page. 

Al.2 Notes concerning chapter 4. 

In transforming from system S' ' to the matrix of 

transformation was ( SBA ) and the position of the moon in the 

SE system was found to be 

~r = SBA “r • ’ Al. 3 

se ■ 

The vector 'r ' ' has only an x-component in system S ' ’ 
since the moon is on the x' '-axis at a distance r from the 
origin. That is, 
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or represented as a column matrix. 
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Figure Al.l Euler angles 6 , ^ and . The coordi- 

nate transformations between the primed and the unprimed 
systems can be conveniently expressed in terras of the 
Euler angles. (Goldstein, 1950) 


2 



Figure Al.l 


Figure Al.l 
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This column matrix, when multiplied by ( SB&) from the left 
gives 



Al.7 


Hence equations 4-8 . 
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Appendix 2 

nf^riva-bion of eauations 3- 2 through 


Gonsider a charged particle inoying in a region in which 
there exists electric and magnetic fields . The force on the 

particle due to the electric field is, -v:; 


F = q E 
e ^ 


A2.1 


and the force due to the magnetic field is 


F, = q V X B 
b 


so that the equation of motion (eq. 3-1) is 

dv 


A2 .2 


m 


dt 


or in component form , 


dV 

^ dt^ 


dv 

^ dt^ 


dV 

dt 


"e is always perpendicular 


= q + 

“v x 'b ) 

A2.3 

! + V B 

X y 2 

- V B ) 
z y 

A2.4 

I + V B 
y z X 

- V B ) 

A2.5 

I + V B 

- V B ) 

A2.6 

z X Y 



given by 
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WB can taks B=Bz and. E~Ey. Equations 4f5/ and 6 of this- 
section then become, 

m -^x = q V B A2 . 8 

dt y 

III = q E - q V B A2 , 9 

dt-^ X 

m-|rz = 0 A2.10 

dt 

Taking the time derivative of equation A2.8 and substituting 
it into equation A2 . 9 gives 

2 

m(m/qB) -^x = q E - q B A2.11 

or 

2 

-^x + (qB/m)^V =(q^ E B/ra^> A2.12 

dt^ ^ 

Defining 

tii = q B / m 

A2.13 . 

and 

. 2 "..■■■■ ■ 2 

K = (q E B/ ra ) 

equation A2 . 12 becomes 

■'.'v 2 

... j \7 O' - ■ ■ 

.^-^x + V = K A2.14 

at^ 
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Equation A2.14 is a well known second order, first degree 
linear differential equation with constant coefficients whose 
solution is of the form 


V (t) = A cos wt + D sin &Jt + c 


A2 .15 


The cons tant G can be determined by subs titutihg equation A2 . 15 
into equation A2 .14 . i .e . 


V = 

X 

V = 

y 


- A sin <ut + D cos (lit. 

2 2 

- Acoscjt - Dsin (ii t 

2 

- 6> (V - C) 

X 


A2.16 


SO that 


becomes 


or 


V + CO V = K 

X X 


C CO = K 


A2.17 


A2.18 


C = i^co = E / B 


A2.19 


A and D can be determined by applying the boundary conditions 
at t=0 and equation A2.8. 


A = V - E / B 
xo 


A2.20 


and from equation A2.16 and A2.8 
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(0) = 6) D = qB/m V 

X yO 


A2.21 


or 


D = V 


yo 


A2.22 


and hence. 


V (t) 

X 


(V „ - E/B) coscjt + V sin (Ut + E/B 

XO yo ' 


Similarly, 


Vy(t) 


cos Cjt " “ E/b) sin cot 


A2.23 


A2.24 


arid from eq . 2 Jl2 . 10 , 


( t) = V 
z zo 


A2.25 


The quantity E/B will be shown to be the magnitude of 
the drift velocity. Therefore the equations for the velocity 


are 




= (V - V-) cos cot + V sin cot + v, 

xo - - d yo . v;:.. ' , d 

= V coscot - (V - V^j sin CO t 

-.yo./,-.,.;:. ■ xO"-, ■:■■■, d 

'I... 

.■ , ZO . 


A2.26 

A2.27 

A2.28 


Integrating eqs . 26 through 28 of ’ this section arid 
plying the condition that the particle start cut at the ori- 
gin yields. 
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X(t) 

y(t) 

Z(t) 

Equations 29 and 30 are the parametric equations of a cycloid 
in the x-y plane and the z~motion is independent of the x-y 
motion. The transformation from a fixed coordinate system to 
one that is moving parallel to the x-axit, v?ith velocity Vpis 

X - V^ t 
d 

Y A2.32 

Z 


X' = 
Y' = 
Z‘ = 


= j(V -V,) /fu j sin w t + (V /<u ) (1-cos £i) t) + V.t A2.29 


xo d 


yo 


= (V /cu) sinwt + /(V -V.) /<u) (cos tu t - 1) A2.30 

yo V xo d / 


= V t 
zo 


A2 . 31 


This implies that eqs. A2.30 and A2.31 are unchanged in 
such a transformation but equation A2.29 becomes , 

X (t) = /(V^^—V^) /<u ^ sin <ut + V^^/ta (1—Gos tut) 


Let A = (V -V )/o) and d = V /tu 
xo d yo 


9 !>=(Ut then. 


X' = A sin 9 ^ + D (1 - cos 9 !) ) 
Y ' = D sin 9 !) + A (cos^ -1) 


A2.33 


Rewriting eqsi A2,33, 
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(X' ~ D) - A s±n<p - D cos<p 
(Y' -~("A)) = D sinijJ) + A cos<p 


A2.34 


Defining the angle r 


6 = tan"^ (A/d) 


A2.35 


then equations A2. 34 become. 


' D) = "^/a^ + D^ cos ( 6 + i> ) 



(Y' -(-A)) =\Ja^ ri- sin( $ ) 


A2.36 



Let R =^/A^ + , then 


(X ' - D) ^ + (y - (~A) ) ^ 


A2.37 


Apparently the equations A2. 29 and A2.30 represent cir- 
cular motion when transformed to a system that is moving par- 
allel to the x-axis at velocity V^. The center of the circle 
in the moving frame is , 


X 

c 


D = V /o) 
yo 



(V, 




ft2,38 



and it has radius/ 
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R = l/co ( 


(V -V,) 
xo d 


+ V 


, 1/2 ^ . 

) ' =m^ / qB 


yo 


A2 .39 


where is the perpendicular velocity in the moving frame. 
Hence equations 3-5 and 3-6 . 

The particle appears to move under the influence of a . 
magnetic field only in the moving frame. That is, we have 
transformed into the frame in which there is no electriC field. 

= y{ + V X B ) A2.40 

Consider the relation, 

^ X (%X 'b) = B A2.41 

X "b) 23/B^ A2„42 

but from eq. A2.40, 

- ~V X "B A2.43 

(EX^)/B^ ' 

V7hich is equivalent to equation 3-8. 

The magnitude of V is ■ ■ . ' 


Evidently) 




V = 


E, / B = V- A2.45 
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lienee the particle obeying eg. A2. 3 seeks circular motion 
in a frame in which the electric field is sero, i.e. one 
moving with velocity "v such that the component of the vel- 
ocity that is perpendicular to the field is equal to the 
drift velocity. Obviously, the particle can move in the z- 
direction as well as in the x~y plane, but any velocity it 
has in the z-direction is velocity parallel to the field and 
it will still be moving in a frame in which the electric 
field is zero. Again, since the z-motion and the x-y planar 
motion are uncoupled, the total motion can be described as 
cycloidal in the X-y plane and uniform translation of that 
motion in the z-direotion. The resultant velocity vector of 
the particle's guiding center is 

*V = V z + X A2.46 

r zo d 

which is equivalent to eq. 3-7. Clearly, the particle moving 
in a frame with velocity"^ sees no electric field. 

■^ = y ( *^+ (V^^z-i- x) X B z) 

= y ( V^B X X z) A2.47 

Si-Xice the perpendicular E field is due to, the solar wiud 
velocity in the present treatment, using eg. A2.44 and 
B it is easy to see that 
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V = V - (V » B) B A2.48 

d SW sw 


= ( -(V X B) X ^ )/ B^ A2.49 

d sw ' 


or using the identity, 

A xCb X c)= (B»A) "C - (C-A) ^ A2.50 

■y = (V B^ - (V * "B) ■! ) / B^ 
d SV7 sw 

V -V - {V •%) B A2.51 

d sw s'v; 


hence equation 3-9 • 



Appendix 3 


a.3 « 1 A representative correlation prograiri is discussed in 
this section, viz. the program that produced figure 5-6. The 
basic idea behind all of the correlation programs was discus- 
sed in chapter 5 and will not be repeated here. Only the de- 
tails of the procedure used to obtain the results are relevant 

to this discussion. 

Program A3.1 does the follov;ing; 

a. ) reads the Explorer 35 magnetometer tape to find the 
start time 

b. ) reads the SIDE data tape to find start time 

c. ) reads a record of Explorer 35 data to determine the 
value of the field latitude for that 82 second interval and 
decides the appropriate bin number for that latitude (i.e. for 
the 180 values of field latitude from -90" to +90", 90 bins 

of 2" each are used and are numbered from 1 to 90) . 

d. ) after selecting the bin number, the average integral 
flux over the same 82 second interval is calculated from the 
SIDE data tape and this number is deposited in memory at that 

bin location . 

e. ) The processes c. and d. above are repeated until the 
period of interest has been evaluated. Each time flux is ad- 
ded to a particular bin, a counter for that bin is incremented 
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so that the number of times that bin was enGoimtered is 
known . 

f .) at the' end of the period of interest, the total 
accumulated flux in each bin is divided by the number of 
times that that bin was encountered. Thus the number ob- 
tained in that bin represents the average integral flux 
observed at a given latitude over the period of interest. 

g.) The results of all 90 bins are then dumped onto 
magnetic tape and are subsequently plotted to obtain figure 
5-6. 


The details of the calculation of the integral flux are 
as follows. Figure A3.1 shows the meaning of the parameters 
relevant to the following discussion and the approximate na- 
ture of the integration over the energy spectrum. The appro- 
priate values of the parameters for the Apollo 14 SIDE are 
given in table A3.1. 

For the ith energy channel, the number of counts C- will 

be given by , 

e = R. At R. = count rate (sec 

i i 2. 

At = 1.13 sec 

the count rate R. is related to the differential energy flux 


by 


Figure A3.1 Illustration of the technique used to cal- 
culate integral flux. The nuniber that represents the 
integral flux is given by summing the areas of the unshad 
ed rectangles defined by AE^ and (dj/dE) . 



Figure A3 



Table A3.1 


nominal energy 

peak response 

AE 

•5 2 

cm “ 

G 

str) 

Ae (fwhm) 

1 

10 

7 

8.5 

3.6 

1.4 

2 

20 

17 

11.5 

5.0 

1.0 

3 

30 

30 

14.0 

6.8 

2.0 

4 

50 

45 

19.5 

3.2 

3.0 

5 

70 

69 

25.5 

15.0 

6.0 

6 

100 

96 

88.0 

13.0 

6.0 

7 

250 

245 

214.5 

7.7 

24.0 

8 

500 

525 

242.5 

6.5 

30.0 

9 

750 

730 

237.5 

3.3 

30.0 

10 

1000 

1000 

235.0 

9.2 

102.0 

11 

1250 

1200 

225.0 

4.5 

50.0 

12 

1500 

1450 

262.5 

3.2 

76.0 

13 

1750 

1725 

287.5 

9.0 

140.0 

14 

2000 

2025 

262.5 

4.7 

230.0 

15 

2250 

2250 

237.5 

7.9 

180.0 

16 

2500 

2500 

225.0 

9.8 

190,0 

17 

2750 

2700 

225.0 

3.3 

250.0 

18 

3000 

2950 

262.5 

9.9 

320.0 

19 

3250 

3225 

287.5 

8.8 

330.0 

20 

3500 

3525 

262.5 

7.1 

300.0 



154 



where , 2 

Gi = G^(E) (cm -str) 

Ae^ = energy passband (PWHM) 

The flux, integrated over the energy spectrum is j. 



= width of the i th integration interval 


or in terms of the known parameters, 

20 / A 

J I So 


J = \ ^ "i _ 1/2 

At Ae . G . 
i=l 1 1 


At Ae^g^ 


At Ae g 
20 20 


(See figure A3.1 to understand the second term) 

The values of were determined by taking half of the 

energy difference between the i th peak and the i+1 th peak 

and half of the difference between E. and E, . , i.e., 

1 1— J. 


>i = { ^i+i - =i-i} 
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A3 .2 Program A3. 2 is the program used to generate the syn- 
thetic data of chapter 5. The theory underlying the gener- r 
ation of these data was covered in chapter 4. The program 
presented here is the SDS 910 version. The IBM version is 
somewhat lengthy and will not be included. It is basically 
the same as the 910 version except that the ephemeris data 
are calculated in that program. However, the tests relevant 
to the acceptance criteria are the same. 

To facilitate' the understanding of program A3. 2, it should 
suffice to define the parameters relevant to the theory of 
chapter 4 . 

T=kinetic energy of the particle 

D=d=distance from the center of the earth to the nose 
of the shock . 

CAPD=D=distance from the center of the earth to the wing 
of the shock measured along the y-axis. 

-CR=r=position vector of the moon 

APX,APY,2iPZ,=components of the detector look vector in 
a solar ecliptic system d 

VPx=V =resultant velocity vector for the guiding center 

X 

VDOTAP= V • b 

X 

B2 and AC4 represent equations 4-22 

HGT=counts 


C IRY< 2S) « BIN t 91 J» IT.U£( 5) . IDI V( 91 ) .FLX ( 10 ) . IWt 20) 

CALL ITYPE(ITIHE) 

C TYPE I M START TIME 

CALL I TY PEC I RUM) 

C TYPE IN RUNTIME IN HOURS 
CALL ITYPEC IW) 

C TYPE IN HODIFIEO CEOHETRIC FACTORS 
HT0DAY=1./24. 

T0DAY^HTaDAYf6O. 

ST0DAY=HT0DAY/36\, , . 


1 


OGAP=i .fTGDAY 

TB=FL0 ATUTT 2 r ) +RLC ATI t ITI-HE ( 3 ) ) »FTDDA Y+ FLDA T( IT IM G< 4 > ) *TCDA Y+FL 
IQATC IT IHEC5) r»STQOAY 
TEND=T8+RUNT IK 


UD 1 K = l,9 1 
aiN(KJ=0. 


IGIV<K) = : 

ASSIGN 2 TO IPAR 


ASSIGN 3 TO IE OF 


I CONTINUE 

CALL SETUPC 1. 28. IRY, IRECU IPAR, lEOF ) 

C FIND START TII-E ON EiXPLORER JS TiAFE 
17 READ INPUT TAPE 2 , 16 , I D A. Y , I H R , If' I N , ISEC 
16 F0RHAT(4I3T 


18 

19 

20 
21 


IDAY=I0AY+1 

IFCIDAY-ITIHEI 2) M 7. 18.2 i 
IF( I rfl-lTIHEl J >) 17. ) 9. 21 
IFC IMI N- IT IHEC4) ) 17. 2 C , 2 I 
IF ( ISEC - IT IHEC 5) ) 1 7. 21.2 . 

TB=FLO AT( 10AY)i-FLGAT( I HR MHT.OEAY LCATC IMIN) »T OC A Y+FL OAT U SEC ) »S TL 


IDAY 

C FIND START TIME ON SIDE TAPE 
SI CALL INPUT 

II T=FLOAT( IRYC 2) ) + FLDA T 1 1 R Y C 3 ) > ♦ HTDC A Y+FL OA T ( IR Y I 4 I ) *T OD AY i-( FL OAT I IR 

lYl 5 > »/ JU. ) fSTOOAY 


25 IFCT-TB) 23.2 4, 24 

2 3 T=T+2 4. 16»ST0DAY 

CALL INPUT 
IFC IRYC 2 ) ) 25 ,22, 22 
24 TPLUS=T + 24 .16*ST0DAY 

TLAST'TB 

C 8CTH TAPES ARE ON START TIME HERE 
C READ EXP 35 TAPE ^ 

r READ INPUT TAPE 2. 2 6 , 1 0 A Y . IHR , IM I N , ISEC, BL AT 




OF POOR QUALffi 


A3.1 
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26 


IFCRHAT(*Ii*32XtF6.U 

TB =FLO Airi DAY) + F LC A T C I HR ) »HT;a C AY + F L 0 A T < IH I N > 


»T ODA Y+ FL OAT 1 1 S£ C ) *S T C 


lOAY 

gap=tb-tlast 
IF(CAP »30, 30,5') 

5T- IF{6AP-0GAP) 48, A 8, 51 

C CALCULATE BIN NUMBER 
48 I=BLAT/2. 

1=1+46 

N=0 

SUH=0. 

C CALCULATE INTEGRAL FLUX 
29 CALL INPUT 

N=N+1 
CT=0. 

DO 4 J = 1»Z0 

IF{ IR Y t J+5 > ) 6 .,4 ,4 


6 ' 

4 


CT = CT+FLDAT{ IRY{ J+5) J*FL DOT ( IW( 21- J) 
CT=CT-{ FLOATl IRY t6 > ) » I Mt 2 - > + FL DA T( IR Y< 25 ) ) 


*IH( 1 ) )/ 20- 


62 SUM=SUH+CT 

IF ( IRY { 2) >27.28, 28 


27 TP=TPLUS 

ea ??,FL0«(ll!Y(ai..FLDAT(18Y(3n.HTC0«.FL0ATUR»<4)l.T001Y*<1=L0»t(I 

IRY (5 ) > / I . - 0. » <STQDAY 

33 TPLUS=TP+24. 16»STDDA Y 

IF(TP-TB >29 . 11, 3 I 

31 SUM=SUH/FLOAT (N) 

BIN! I»=BIN(l J+SUM 
lOIVU )=lDIVtl)+l 
TLAST=TB 

[F (TB-TEND) 30, 32 .32 

32 CONTINUE 

CALL I TYPE! I MATE » 

C MCUNT BLANK TAPE ON UNIT 1 
iDC 13 K= 1, 81 . 10 
DO 5 7 J = 

I=K+ J- 1 

IF { I'D IV< I J > 14, 14, 56 


14 FLX<J>=C. 

GO TO 57 

56 FLX{J)=BIN(I )/FLQAT( ICIVKI)) 

57 ICCNTINUE 

^ ““"s«rT“ouTruT ■ 

i\KiirNunUriiF°rlMb ?’ aa"c luv”')'* 

13' WRITE OUTPUT TAPE 1,58, ( ICIV t- ), J- K,K V) 

58 format <8X, 10 I l'' l 
END FILE 1 

3 ICCNTINUE 

END 


>g ®3 
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C PROGRAM A3. 7 

C J. BENSON RICE UNIVERSITY. 

C IP(U = DIST ANCE TO NOSE OF SHOCK UN RE 

C IP(2 )/ IP< 3 ) = DIST ANCE TO MOON- AT BOW SHOCK CROSSING IN RE 
C IF(4)/HFt3) = SCLAR ECL LONG AT BOW SHOCK CROSSING IN RADIANS 
C IIF(6)=S0LAR WIND VELOCITY IK» I'M / SEC 
C IP(7)=SI2E OF THE LOOK CONE IIA DEGREES 
C IP(8)=SHA00W ANGLE I N DEGREES 

C ;SFACOW ANGLE IS USED. TO lAVOID UNECESSARY TESTING 
C IP(9 )=NUMBER OF RECORDS iTC BE PROCESSED AT A T IME 
C IFd )=SOURCE NUMBER 

DIME NS ION T( a ; IP( lu) , IRI 24 I ,R< 12 ),1CUT ( 28) 

EQUI VALENCE! IR ( 2 J .RC i ) ) 

ASSIGN 40 TO IPAR ‘ 

ASSIGN 4 I TO I EO F 

C CALCULATE 2‘ VALUES OF THE ENERGY 
IDC 2 K = 1.I4 

2 T( 21-K )= 3500 (K-i )» 250- 
T< 6) = I' T. 

IOC 3 K = 1.3 

3 r(6-KJ=7 .-(K-l)*2^. 

T( 1) =1'". 

T(2>'=2T.. ■ . 

CALL ITYPE(IP) 

C TYPE IN INPUT PARAMETERS 
0- FLOAT! IP! 1 )) 

RE = FLnAl! IP! 2)1/ FLOAT! IP ! ;1 ) 

SELAHS=FLQAT !IP( 4 )) / FL OA T ! IP ! 5 ) ) 

VS=FLOATC IP! 0) ) 

SLA=FLOAT! IP! 7 V) 

SHAO=FLOAT! IP!8 ) ) 

BETA=. 8726 
PI =3. 14159 
PI2=PI n. 

ITCRAD=PI/180- 
SL A= SL A* TOR A D 
CLA = SIN! PI2- SLA) 

SeS-SlN! SELA8S) 

CES = SIN! P12-SELABS) 

C CALCULATE VALUE OF H IWG D I ST Af.CE C APO 

CAPD=SQR T! C-O'^RE’REv saS^.SBS ) /! RE^fC 8S-0 ) ) 

SB=S IN ( BETA ) 

CB = SIN !P I2-B ETA) 

A" CALL SETUP! 1 .2^. IR. I REC. IPAR , lEOF ) 

IOC 7 i = 1 , IP! 9) 

C READ EPHbMERIS DATA 
call I NPUT 
7 CONTINUE 
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C TEST SHADOW ANGLE 

IF(R { i n-SHADM, 1. 5 

i ■ DG 6 J=l,20 

6 IDUT(J+5)=0 

■GC TU 4 

S IF (R U 1 ) -1 80 . )9, ifa, 1 6 

16 IF(R(i U-SHAO-18,;. 11.1.9 
9 ICCNTINUE 

C CALCULATE SOLAR ECLIPTIC ANGULARl PARAMETERS 
SELD=S IN<R(7 ) ) 

CELO=S INC FI2-R( 7 ) » 
ssELDD=s INCR tan 

CS ELOD = S IN(P I2-R C8 ) ) 

SELA=S INCRC 9 )) 

CELA=S IN(PI2-R(9 J J 
SSELOA=SIMR (10) ) 

CSELOA = SIN(P I2-R ( 1^) > 

S8L0N = SIN(R( ll)*TORAD) 
t6LDN= SI N( PI 2-RC ID'-TCRAO) 

SOLAI^SINCRC .2)»T0RA0) 

CBLAT = SIN(Pl 2-R( i2)»TCRADI 

BX=CaL AT »CBLOW 

BY = CBL AT-fSBL UN 

8Z = SBL AT 

IF (BY I 17, 19, 19 

19 BX=-8X 
8Y=-8Y 
BZ =-BZ 

C CALCULATE R 

17 CR = S QR T C R( 4 ) »R( 4 ) +R ( ;> ) *R ( :) ■'•R < B S *R ( ^ * 
CRX = CR»CELA*CSEL UA 
CR Y = CR »CEL A» SSEL UA 
CRZ=CK »SELA 

C calculate DETECTOR LUCK \/ECT>OR 
AP X=-C EL D»CS EL 00 
APY = -CELD»SS6L0D 
APZs-SELD 

DUNG = VS'‘CE'‘CBL AT»CBU0A-VS*SB »CRL AT *SBL0N 

C EVALUATE 2 ENER GIES 
DO 4 K = 1 , <:0 

iVF = S(lRT< Tt K ) * 1.6 2 E - 22 ) » i. "9 36 El 2 
HGT= .. 

C EVALUATE i ^ PITCH ANGLES 
DO 13 J= 1 , 10 
PA = ( J- 1 » i 0. » TOR AO 
bP A=SIN( PA ) 

CPA=SGRTU SPA» SPA) 

VP =VP* CP A 


Pmgmm; page is 

lOl.lPOOB ([QUALITY 
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C CALCULATE RESULTANT V6LDC1TYI \iECTCR 
VRX = -VS»CB + BX»OUNG<-BX*VP 
VRY=VS*SB+BY»UUN G+BY »VP 


VRZ=B2*(DLNG+VP) 

VR=SQRT( VRX* VRX+ V«Y* VR Y+ VR2* VBZ) . 

C CALCULATE ANGLE BETWEEN VR ANC OET LOOK VEC AND TEST 
VDQTAP=< VRXf APX* VRY* APY+VRZ* APZ) / VR 
IF (VDGTAP-CLA> 13,10, K 
in IF(VDQTAP-1. )18f i3 . 1 3 
C TEST FOR SOURCE POINT 

1/VRY n 

PARF=I(CRY*VRZ>/VRY» »( (CRY»VRZ)/VR Y) 

BELCH=-CRY»( ( ZoCR2*VR2 > /IVHY+ t C APO* C APO* VRX 1/ (D^liR Y) ) 
6RUCT={C AP0» CAPO ■'‘CRX ) / D+ CRZ» CRZ-CA PO *C APO 
CAPC = B AR F+BELCH+ERUC T 


Bi=C APB»CAPB 
AC4 = ‘,»CAPA»CAPC 
IF (82- AC^) 13 , 1 1, 1.1 

n IF (PA) 13, 14, 15 

C FILL CHANNELS WITH APFROPRIATt COUNTS 
U H(jT-IP( l J)*l .♦TORAD 

GO TO 13 

15 hGT= HGT+ IP( 1 1- 1 *5 IN t PA) * 1 ■»TQR AC 

IJ ICENTIAUE 

I0UT(2 6-K)=HGT 
4 CONTINUE 


IOUT( 1 )=-l 
'1CUT(2)=R( 1> 

IOUT< 3 )^R( 2) 
lOUT (4 )=R( 3) 

!1CUT(5) = J 
iDC 8 J = l.i 

^ CALL SETUP! 2 . 28.1CUT . IRE C, IP AR , lEO F ) 

C WRITE A SPECTRUM ON MAG TAPE 
CALL OUTPUT 
GO TO 4r 

41 CONTINUE 
END 
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Appendix 4 


Brief review of the acceleration mechanisms. 

A4 . 1 Sonnerup model - The basic idea of the Sonnerup 
model is simple specular reflection from the shock and sub- 
sequent energization of the particle by the interplanetary 
electric field. To get an order of magnitude estimate of the 
increase in energy after reflection we can calculate the ra- 
tio of the reflected energy to the initial energy. Assuming 
the reflected energy is given by the initial energy plus the 
energy given to the particle by the electric field over one 
half of a cyclotron orbit, 


f = 1/2 mV + q n r 

r sw ^ IP 


A4.1 


where 


r -■= mVg^qB 


A4.2 


and 


^ = angle between V and B 

sw 

2 2 

f = 1/2 mV + n-mV sin^ A4.4 

r sw sw 


so that. 


f /f. = 1 + (ravj 77sin0)/(l/2 mV^ ) 
r 1 sw sw 


A4.5 
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f /e. = 1 + 27TsiniA A4.6 

^ X 

for a typical value of ^ , ( ^ = V4) > 

fr/fi == 5-45 A4.7 


That is, a particle could gain roughly 5.5 times its incoming 
energy by this process. For a typical bulk flow velocity of 
about 400 km/sec, the reflected particle energy would be, by 
eq. A4.7, about 4.7 keV. This is the right order of magnitude 
for the observed particles in the upstream region. For dif- 
ferent values of ^ the energy gained varies from 0 to about 
7.28 times the initial energy. 

Sonnerup treats this problem in a more rigorous fashion 
in his 1969 paper where he considers that the particle energy 
is conserved in the frame of reference in which the electric 
field vanishes . The results of this more rigorous treatment 
will be stated here. By considering the reflection process in 
the frame of zero F field and the rest frame of the shock, 
Sonnerup obtains. 


/ _ , n (1+5) cos { 0+ y ) sin(^ - (jf*) sin{ + y ) _ ^ 

f/f, — X + Z - r A4 . 8 

r X ( 1+ /X ) 003“=^ ( 0 - 0 ) 

where the angles are those shown in figure A4.1 and the def- 
initions of IX and S are as follows. 


2 , 2 
V_L./V, 


S= V,, / V,,. 
'•r Hx 


A4.9 
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Figure A4.1a . Vector diagram showing how the incident 
particle guiding center velocity vector, , xs decom- 
posed into a component, *V,j ^ , along the field and a com- 
ponent along the shock front. 

Figure A4.1b. Similarly the velocity vector of the re- 
flected particle is decomposed into ^ along the field 
and along the shock front. The angle ^ betv?een the 
field direction and the solar wind vector is the 

garden hose angle. The angle between the shock normal 

and "v measures the obliqueness of the shock at the point 
sw 

of reflection. The angle V is non zero only if the in- 
cident particle has a non-zero thermal velocity component 
parallel to the field. 

Figure A4.1c. Ratio of the particle energy after and be- 
fore reflection as a function of shock obliquity. 

[the above is paraphrased from Sonnerup (1969)] 
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Figure A4.lci shows the energy ratio as a function of shock 

obliquity ( 4 ^) for the values (f = y=o /*=o. 

It should be noted that the singularity in equation A4.8 at 

f Vl. corresponds to the case in which transformation to 

the frame in which the E field is zero is not possible. That 

is, the singularity is artificial, physically, and simply 

means that no finite velocity vector “^^(see fig. A4.1) can 

t 

be found that makes the solar wind velocity parallel to "B in 
the moving frame. It can be seen from examining the curve in 
figure A4.I that the Sonnerup model predicts very high values 
of energization on the dusk side of the magnetosphere and rel- 
atively small ones on the dav/n side. Although electrons of 
energy greater than 40 keV have been observed [ Anderson, (1968) ] 
and more recently protons of energy 30 to 100 keV [Lin et al., 
(1974)] in the bow shock regions, no such asymmetry has been 
observed. The dawn side energy increases from the model are 
much too low to account for the observed high energy particles 
in this region. On these bases one must conclude that the 
Sonnerup model is probably not correct. 
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A4.2 Fermi acceleration mechanism- [Jokipii , (1966) ] 

The basic idea of the Jokipii model of Fermi accel- 
eration at the bow shock is that quasi -thermal ized parti- 
cles escape from the shock region into a region that is 
bounded on one end by the bow shock and on the other by a 
magnetic irregularity that forms a mirror point upstream 
of the bow shock, (see figure A4.2) Under the assumption 
that the particle velocity parallel to the field is much 
greater than the velocity of the moving irregularity, the 
number of collisions per second that the particle makes with 
the mirror point is 


V / 2 L(t) A4.10 

P// 


where 

L(t) = - Vt A4.11 

and the gain in momentum on each collision is 

Ap = m V A4.12 

This treatment differs slightly from Jokipii 's but the 
differences are primarily notational; for instance in fig- 
gure A4.2, and t^=0 . 
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Figure A4.2 Jokipii model of the Fermi mechanism for 
accelerating particles at the bow shock. The trapping 
region is defined by the bow shock at one end and a mag- 
netic mirror at the other. The magnetic irregularity moves 
tovjard the shock with velocity 


is 



Shock 

Front 


Figure A4 . 2 
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From equations A4.10 and A4.12, the time 
gain of momentum is 



or 


^ = mV V/2L /2L 


dt p, 


where 


and 


II 


IX = V /V = cos (pitch angle) 

P P 
^11 


but since. 


p = m V 


, 24 ^ 2 2 . 1/2 

W = + P c ) 


2,2 22 

p=W/c 


dp 

dt 


2 W/ c" 


dW 

dt 


rate of 


A4.13 


A4.14 


A4.15 


A4.16 


A4.17 


A4.18 


A4.19 


168 


Comparing equations A4.19 and A4.14 we' see that 

A4.20 

which differs from Jokipii’s result only in notation. This 

differential equation was derived under the condition that 

V < < V . If the magnetic moment of the particle is con- 
P 

served in collisions with the shock and mirror then as the 
particle gains energy its pitch angle decreases and it will 
eventually be lost to the solar wind. This sets a limit on 
the energy that can be gained in a process like this and that 
limit depends on the dimensions of the trapping region and 
the velocity of the irregularity in the field. One gathers 
from this simplistic approach that a particle could gain par- 
allel energy only up to about that associated with the bulk 
flow velocity of the solar wind. Energies observed in the up- 
stream regions are much higher than this (up to 100 keV protons) 
[Lin et al. ( 1974 )] which requires some other process to pro- 
due these particles. Jokipii depends on pitch angle scatter- 
ing at a fluctuating shock to avoid losing his particles too 
quickly. Lin et al . [ 1974 ] have pointed out that this accel- 
eration mechanism predicts a hardening of the energy spectrum 
toward the v/ings of the shock due to the longer time that a 
trapping region ic in contact with the shock V7hich is not ob- 


dW 

dt 


= (AtV/L(t)) W 


2 4 2 

1 - (m'^c /W ) 
o 


served . 


169 


Acknowledgements 


I wish to thank my thesis advisor, Dr. John W. Freeman, Jr., 
for giving ire the opportunity to work with a unique project, 
the Suprathermal Ion Detector Experiment. Dr. Freeman's pat- 
ience and understandxng have heen xndispensatle xn the fxnal 
realization of this thesis. 

I would also like to thank Dr. H. Kent Hills for countless 
hours he has spent discussing the subject matter of this thesis 
with me. His immense knowledge of computers and his astonishing 
attention to. detail have been invaluable in the sucessful com- 
pletion of this work. 

I have benefitted greatly from my professional association 
and personal friendship with Dr. Richard R. Vondrak who has con- 
tributed much of his time to valuable discussions concernxng 
this thesis. \ 

I would like to thank Tom Rich for his help in computer 
programming and Dave Nystrom and John McGarity for maintaining 
the sanity of the SDS 910 computer during the past three years. 

I have enjoyed knowing and working v/ith the other members 
of the SIDE group, past and present, namely. Dr. Martha Fenner, 
my good friend and ex— office mate, Drs . R.A. Lindeman and 
R.A. Medrano, Mohammed Ibrahim, Dave Hardy, Henry Schneider and 


Dr. Jiirg Meister. 


170 


Thanks also go to Steve Bowling and Mike Geis for many 
valuable discussions . 

Finally, I would like to thank my wife Jackie for unself- 
ishly spending many nights as a computer widow and for 
typing the final copy of this thesis. 

This research has been supported in part by the 
National Aeronautics and Space Administration under 


contract NAS 9-5911. 


Bibliography 


171 


Anderson, K. A., Energetic Electrons of Terrestrial Origin 
Upstream in the Solar Wind, JGR , 7^, 2387, 1968. 

Asbridge, J. R., Bame, S. J., Strong, I. B., Outward Flow 

of Protons from the Earth's Bow Shock, JGR_, 73 , 5777, 19 68. 

Axford, W.I., The Interaction Between the Solar Wind and the 
Earth's Magnetosphere, JGR, 67 , 3791, 1962. 

Bowling, S. B., The Magnetic Field Structure of the Plasma 
Sheet at 30 Re, M.S. Thesis, Rice U., 1973. 

Colburn, D.S., Currie, R.G., Mihalov, J.D., Sonett, C.P., 

Diamagnetic Solar Wind Cavity Discovered Behind the Moon, 
Science , 158 , 1040, 1967. 

Colburn, D.S., Mihalov, J.D., Sonett, C.P., Waid, S.H., The 
Lunar Cavity in the Solar Wind, Trans . Amer . Geophys . 

Union , 49 , 234, 1968. 

Colburn, D.S., Description of Ames Magnetometer Data from 

Explorer 33 and Explorer 35 deposited in the Data Bank, 
NASA/Ames Research Center preprint, 1969. 

Dyal, P., Parkin, C.W., Colburn, D.S., Schubert, G., Lunar 
Surface Magnetometer Experiment, Apollo 16 Preliminary 
Science Report, 11-1, 1972. 

Dyal, P., Parkin, C.W., Sonett, C.P., DuBois, R.L., Simmons, G., 
Lunar Portable Magnetometer Experim^ent, Apollo 16 Pre-: 
liminary Science Report , 12-1, 1972. 

Fairfield, D.H., Bow Shock Associated Waves Observed in the 
Far Upstream Interplanetary Mediiain.- JGR , JH, 14 , 3541, 

1969. 

Fenner, M.A. , Magnetosheath Plasma at 60 Re, M.S. Thesis, 

Rice U . , 1971 . 

Fenner, M.A., Freeman, J.W., Jr., Hills, H.Ky The Electric 

Potential of the Lunar Surface, Proceedi ncs of the Fourth 
Lunar Science Conference , 3 , 2877, 1973. 

Fenner, M, A., Observations of Magnetosheath Plasma at the 
Lunar Orbit, Ph.D. Thesis, Rice U., 1974. 

Frank, L.A. and Shope, W.L., A Cinem.ato graphic Display of 

Observations of Lov; Energy Proton and Electron Spectra 
in the Terrestrial Magnetosphere and Magnetosheath and 
in the Interplanetary Medium, Trans. Am. G. U. , 49 , 1 , 

279, 1968. 

Prank, L.A., On the presence of Low Energy Protons (5^E^50 keV) 
in the interplanetary medium, JGR , 75 , 707, 1970. 


172 


Freeman, J.W., Jr., The ALSEP Lunar Ionosphere Detector, Rice U 
preprint, 1969. 

Freeman, J.W., Jr., Balsiger, H., Hills, H.K., Preliminary 
Results from the Lunar Ionosphere Detector, Apollo 12 
Preliminary Science Report , 83, 1970. 

Freeman, J.W., Jr., Energetic Ion Bursts on the Nightside of 
the Moon, JGR, _77, _1, 239, 1972. 

Geiss, J., Beuhler, F., Cerutti, H., Eberhardt, P., Filleux, Ch 
Solar Wind Composition Experiment, Apollo 16 Preliminary 
Science Report , 14-1, 1972. 

Goldstein, H., Classical Mechanics , Addison Wesley Pub. Co., 
Reading, Mass., 1950. 

Greenstadt, E.W., Green, I.M., Inouye, G.T., Correlated Magneti 
Field and Plasma Observations of the Earth's Bow Shock, 
JGR, 73, 1 , 51, 1968. 

Heppner, J.P., Sugiura, M., Skillman, T.L., Ledley, B.G., 

Campbell, M. , OGO-A Magnetic Field Observations, JGR , 72 , 
21, 5417, 1967. 

Hills, H.K., ALSEP/SIDE Ephemeris Computations, Rice U. pre- 
print, 1969. 

Hills, H.K., Meister, J., Vondrak, R.R., Freeman, J.W., Jr., 

The Suprathermal Ion Detector Experiment, Apollo 15 
Preliminary Science Report , 12-1, 1971. 

Hundhausen, A.J., Composition and Dynamics of the Solar Wind 
Plasma, Rev, of Geophys. and Space Physics , _4, 729, 

1970. 

Johnson, F.S., Evans, D.E., Caroll, J.M., Cold Cathode Gauge 
Experiment, Apollo 14 Preliminary Science Report , 185, 

1971. 

Jokipii, J.R., A Model of Fermi Acceleration at Shock Fronts 
with an Application to the Earth's Bow Shook, Astrophvs . 

J. , 961, 1966. 

Jokipii, J.R., Correlation of ;t?'30 keV Electron Pulses and 

Magnetic Fields in the Magnetosheath and Beyond, JGR , 73 , 
2, 931, 1968. 

Lin, R.P., Meng, C.I., Anderson, K.A., 30 to 100 keV Protons 

Upstream From the Earth's Bow Shock, JGR , 79 , _4_» 489 , 1974 


173 


Lindeman, R.A., Recurring Ion Events at the Lunar Surface, 

M.S. Thesis, Rice U., 1971. 

Lindeman, R.A., Observations of Ions from the Lunar Atmosphere, 
Ph.D. Thesis, Rice U., 1973. 

Lyon, E.P., Bridge, H.S., Binsack, J.H., Explorer 35 Plasma 

Measurements in the Vincinity of the Moon, JGR , 72 , 6113, 
1967. 

Manka, R.H., Lunar Atmosphere and Ionosphere, Ph.D. Thesis, 

Rice U. , 1972. 

McKenzie, J.F., Hydromagnetic Wave Interaction with the Mag- 
netopause and the Bow Shock, Planet, and Space Sci . , 18 , 

_1, 1-23, 1970. 

Medrano, R.A., Unusual Solar Wind and Solar Proton Events Ob- 
served on the Lunar Surface, Ph.D. Thesis, Rice U., 1973. 

Ness, N.F., Scearce, C.S., Seek, J.B., Initial Results of the 
niP-1 Magnetic Field Experiment, JGR , 69 , 17 , 3531, 1964. 

Ness, N.F., The Earth's Magnetic Tail, JGR , 70 , 2989, 1965. 

Ness, N.F., Behanon, K.W., Cantarano, S.C., Scearce, C.S., 
Observations of the Earth's Magnetic Tail and Neutral 
Sheet at 510,000 Kilometers by Explorer 33, JGR , 72 , 

927, 1967. 

Ness, N.F., The Geomagnetic Tail, Rev, of Geophysics , 7, 1, 2, 

97, 1969. 

O'Brien, B.J., Reasoner, B.L., Charged Particle Lunar Envir- 
onment Experiment, Apollo 14 Preliminary Science Report , 

193, 1971. 

Rich, F.J., Plasma in the Geomagnetic Tail Observed by the 

Charged Particle Lunar Environment Experiment, Ph.D. Thesis, 
Rice U. , 1973 . 

Scarf, F.L. et al.. Direct Correlations of Large-Amplitude 
waves with Suprathermal Protons in the Upstream Solar 
Wind, JGR , 75 , 34 , 7316, 1970. 

Sonnerup', B.U.O., Acceleration of Particles Reflected at a 
Shock Front, JGR , 74 , _5, 1301, 1969. 

Smith, E.J., Holzer, R.E., McLeod, M.G., Russell, C.T., Magnetic 
Noise in the Magnetosheath in the Frequency Range 3-300 Hz, 
JGR, 72, 19, 1967. 


174 


Spreiter, J.R., Suiraners, A.L., Alksene, A.Y,, Hydomagnetic 
Flow around the Magnetosphere, Plan, and Spac. 

2A, 223, 1966. 

Snyder, C.W., Clay D.R., Neugebauer, M., Initial Result from 
AJjSEP Solar Wind Spectrometer (abstract) , Trans. Am. 
Geo. U., EOS, 51, 407, 1970. 


Snyder, C.W., Clay, D.R., Neugebauer, M., Solar Wind Speedom- 
eter Experiment, Apollo 12 Preliminary S c ience Rsp.Q.?_t / 

75, 1970. 


Van Allen, J .A. , 
of Energetic 
Polar Caps, 


Fennell, J.F., Ness, N.F., Asymmetric Access 
Solar Protons to the Earth's North and South 
JGR, 76, 19, 4262, 1971. 


Wolfe, J.H., McKibbin, D.D., Pioneer 6 Observatxod . a 

Steady State Magneto sheath , Plan, and Spac . Sci . , yo-s, 

1968. 


Wolfe, J.H., Intriligator , D.S., 
with the Geomagnetic Field, 
511-596, 1970. 


The Solar Wind Interaction 
Spac. Sci. Rev. , IQ, Ar 


Zhigulev, V.N., Romeshevskii , E.A., Concerning the . 

of Currents Flowing in a Conducting Medium with the Earth s 
Magnetic Field, Soviet Phv. Dokl. , 4j 859, 1959. 


